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TEMPERATURE  SOUNDING  FROM  SATELLITES 


Compiled  by  Sigmund  Fritz 
National  Environmental  Satellite  Service 
National  Oceanic  and  Atmospheric  Administration 
Washington,  D.C. 


ABSTRACT.   Vertical  temperature  sounding  of  the 
atmosphere  can  be  calculated  from  satellite  radi- 
ance measurements.   The  theory  relating  satellite 
radiance  measurements  to  temperatures  in  the 
atmosphere  for  both  cloudless  and  cloudy  condi- 
tions is  presented.   The  accuracy  of  vertical 
temperature  soundings  from  satellites  in  compari- 
son with  radiosonde  measurements  is  judged  on  the 
basis  of  two  years  of  operational  experience.  RMS 
temperature  "errors"  varied  from  1.1K  to  3K  for 
"clear"  conditions,  and  from  1.6K  to  3.5K  for 
conditions  during  which  high  clouds  were  judged 
to  be  present.   Sounding  instruments  already 
flown  on  Nimbus  satellites  and  instruments  to  be 
flown  in  the  near  future  are  reviewed.   Remaining 
problems  which  require  further  work  to  improve  the 
accuracy  of  remote  vertical  temperature  soundings 
are  discussed. 

1 .   INTRODUCTION 

Quantitative  measurements  of  temperature  at  various  pressure  levels  in 
the  atmosphere  are  among  the  fundamental  observations  required  in  present- 
day  weather  forecasting.  Until  recently,  such  free-air  temperature 
measurements  were  obtained  almost  entirely  from  balloons  and  rockets. 
Since  these  were  available  mainly  from  populated  land  areas  and  a  few  ships, 
large  gaps  existed  in  the  observational  network. 

The  advent  of  satellites  offered  the  possibility  of  obtaining  a  truly 
world-wide  distribution  of  temperature  soundings.   A  method  for  doing  that 
was  suggested  by  Kaplan  (1959).   He  showed  that  vertical  temperature  sound- 
ings of  the  atmosphere  could  be  inferred  from  satellite  spectral  radiance 
measurements  in  the  15|jm  band  of  CO2.   On  the  basis  of  Kaplan's  suggestion, 
two  instruments  were  developed  to  sound  the  atmosphere  from  satellites. 
These  were  the  SIRS  (Satellite  Infrared  Spectrometer)  and  the  IRIS 
(Infrared  Interferometer  Spectrometer).   However,  it  was  not  until  these 
instruments  were  launched  on  Nimbus  3  in  April  1969,  that  the  accuracy  of 
vertical  temperature  profiles  retrieved  in  operational  practice  could  be 
determined. 


Operational  experience  with  SIRS  indicates  that  much  can  be  done  to 
retrieve  vertical  temperature  profiles  for  use  in  operational  weather  fore- 
casting.  However,  many  problems  remain. 

The  transmit tances  of  the  atmosphere,  not  only  in  the  CO2  band  but  also  in 
the  water  vapor  and  ozone  bands,  were  brought  into  question.   The  absolute 
value  of  the  calibration  of  the  satellite  instruments  has  some  uncertainty 
amounting  to  perhaps  1  or  2  percent.   The  influence  of  clouds  varies,  rang- 
ing from  annoyance  to  complete  obstacle  depending  on  the  amount,  opacity, 
and  height  of  the  clouds.   It  became  apparent  also  that  the  departure  of  the 
initial  estimate  of  the  temperature  profile,  or  of  the  so-called  "first 
guess",  from  the  "true"  temperatures  affects  the  final  result;  and  this 
influence  depends  somewhat  on  the  method  of  temperature  retrieval. 

This  report  summarizes  operational  experience  in  the  United  States, 
although  occasional  reference  is  made  to  work  in  other  countries.  As  shown 
in  the  list  of  contents,  the  report  begins  with  the  theoretical  and  techni- 
cal basis  for  temperature  soundings  from  satellites.   The  theoretical  basis 
for  temperature  soundings  in  a  cloudless  atmosphere  is  presented  in  Sections 
2  and  3.   The  troublesome  problem  of  removing  the  influence  of  clouds  and  an 
assessment  of  the  role  of  aerosols  is  given  in  Section  4.  The  theoretical 
basis  is  followed  by  a  summary  of  existing  experience. 

Operational  experience  with  SIRS,  including  an  error  analysis,  is  summar- 
ized in  Section  5.   However  this  does  not  contain  summaries  of  all  the 
research  experience  obtained  from  SIRS  and  IRIS,  some  of  which  had  its 
impact  on  the  operational  procedures.   Discussions  about  the  instruments 
used  so  far,  and  a  word  about  the  instruments  which  may  be  used  in  the  near 
future,  are  included  in  Section  6.  Finally,  in  Section  7  some  problems 
which  remain  before  further  improvement  in  the  accuracy  of  temperature 
retrieval  can  be  achieved  are  discussed. 

The  main  part  of  each  section  was  prepared  by  the  following  people  of  the 
National  Environmental  Satellite  Service,  NOAA,  although  there  was  inter- 
action with  many  individuals.   S.  Fritz  was  the  convener,  or  editor. 

Section  1.  -  Introduction  (S.  Fritz) 

Section  2.  -  Theoretical  Background  (D.Q.  Wark) 

Section  3.  -  Inversion  Procedures  (H.E.  Fleming) 

Section  4.  -  Clouds  (W.L.  Smith)  and  Aerosols  (H.  Jacobowitz) 

Section  5.  -  Operational  Experience  (S.  Fritz) 

Section  6.  -  Sounding  Instruments  (D.T.  Hilleary-IR; 

J.C.  Alishouse -Microwave) 
Section  7.  -  Some  Remaining  Problems  (mainly  Fritz,  Fleming,  and 

Smith) 

Additional  suggestions  were  received  by  mail  from  several  scientists; 
some  of  these  suggestions  have  been  incorporated  into  the  text. 


Mailed  suggestions  were  received  from 

M.T.  Chahine,  formerly  of  the  Jet  Propulsion  Laboratory,  Pasadena, 

California,  now  at  American  University  of  Beirut; 
R.  Fraser,  and  V.G.  Kunde,  both  of  Goddard  Space  Flight  Center,  NASA; 
L.L.  Stowe,  Jr.,  University  of  California,  Los  Angeles. 

Useful  comments  were  also  supplied  by 

C.M.  Hayden,  A.W.  Johnson,  D.S.  Johnson,  and  H.W.  Yates,  all  of  NESS, 
who  between  them  reviewed  the  whole  manuscript. 

2.   THEORETICAL  BACKGROUND  FOR  TEMPERATURE  SENSING  THROUGH 

A  CLOUDLESS  ATMOSPHERE 

The  simplest  case  to  consider  is  the  temperature  retrieval  in  a  cloudless 
atmosphere  resting  on  a  blackbody  surface.   But  even  in  this  case  atmos- 
pheric composition  and  transmittance  must  be  known  accurately  for  certain 
methods  of  temperature  retrieval  from  satellite  measured  radiances.  We 
generally  assume  the  C02  is  mixed  well  enough  in  the  atmosphere  so  that  its 
composition  can  be  taken  as  invariant  with  time  and  place;  this  is  obvi- 
ously an  approximation,  but  probably  does  not  introduce  serious  errors  on  a 
global  scale.   However  amounts  of  other  gases,  such  as  water  vapor,  are 
highly  variable  and  modify  the  atmospheric  transmittance,  expecially  through 
the  lower  atmosphere,  even  in  the  CO2  emission  bands.  The  atmosphere  is  a 
mixture  of  gases;  the  transmittance  of  its  varying  combinations  must  be 
approximated. 

The  satellite  instrument  spectrally  measures  the  emission  from  this  com- 
bination of  gases.   However,  no  instrument  has  a  pure,  monochromatic  trans- 
mission, so  the  spectral  characteristics  of  each  instrument  must  be  measured 
in  the  laboratory. 

Since  the  measured  radiances  depend  on  the  air  temperature,  on  the  spectral 
transmittance  of  the  atmosphere,  and  on  the  instrumental  characteristics,  it 
is  necessary  to  consider  the  relationships  between  these  quantities.   From 
these  relationships,  we  may  be  able  to  derive  the  atmospheric  temperature 
profile  from  the  measured  radiances  and  a  knowledge  of  other  quantities  by 
methods  which  are  not  based  entirely  on  empirical  regression. 

A.  General  Considerations 

The  response  of  a  satellite -borne  instrument  will  depend  on  the  radiance 
of  the  earth  within  the  field  of  view.  The  response,  whether  linear  or  not, 
can  be  determined  from  a  calibration  procedure,  using  a  uniform  source  of 
known  radiance.  The  earth's  radiance  is  generally  not  uniform  (mainly  be- 
cause of  variable  cloudiness),  but  the  response  of  the  instrument  will 
nevertheless  be  a  measure  of  the  radiative  flux  at  the  detector.  The  equiv- 
alent radiance  of  the  calibration  source  is  the  quantity  frequently 


referred  to  as  the  earth's  radiance. 

The  physical  basis  for  sounding  the  cloudless  atmosphere  is  discussed  in 
this  section.   The  problems  introduced  by  variable  cloudiness  will  be 
treated  in  Section  4. 

The  observed  "radiance"  in  a  narrow  spectral  interval  centered  at  v  can 
be  expressed  as  J 

tt/2     2fT 
KvJ   -  J         J      [Kv,,cp,co')  R(v4,cd,cp')]  dcpdcp'  +  E(v.),  (1) 


where  I  «  radiance, 

cp, ©'  m   spherical  space  coordinates  (azimuth  and  elevation, 
respectively), 
R  ■  instantaneous  and  normalized  field  response  of  the 

instrument,  and 
E  ■>  error  of  measurement. 

The  radiance  is  an  ever-changing  quantity  as  the  satellite  progresses  and 
scans  from  side  to  side.  Variations  arise  from  clouds,  different  atmos- 
pheric and  surface  temperatures,  random  errors  of  measurement,  and  changes 
in  concentrations  of  gaseous  and  particulate  constituents. 

Consider  an  idealized  situation  in  which  there  are  no  clouds  in  the  field- 
of-view  and  temperatures  are  uniform  over  the  horizontal  field.   Then  the 
measured  radiance  is,  aside  from  errors  of  measurement,  proportional  to 
I(v>cpjCp').   In  this  degenerate  condition  the  equation  of  radiative  transfer 
may  be  approached. 

In  an  idealized  atmosphere  which  is  plane-parallel,  is  free  of  scattering 
agents,  and  is  in  local  thermodynamic  equilibrium, 

dKv.coV)  =  [-I(v,co,cc')  +B[v,T(s)]  j  k(v,s)p(s)ds  ,         (2) 

where  B  =  Planck  radiance, 
s  =  geometrical  path, 
k  =  mass  absorption  coefficient, 
p  =  density  of  absorbing  gas, 
v  =  spectral  frequency,  and 
T  =  temperature. 

By  means  of  the  hydrostatic  equation,  the  variable  "s"  can  be  transformed 
to  pressure  so  that  pds=.  £  sece  dp.  the  path  of  the  beam  is  at  angle  e> 

g 


from  the  local  vertical,  "p"  is  atmospheric  pressure,  "qn  is  the  mass  frac- 
tion of  the  absorbing  gas,  and  "g"  is  the  acceleration  of  gravity.  Eq.  (2) 
now  becomes 

dl(v,e)    =   [+Kv,  9)    -  B[v,T(p)]  j    k(v,p)   iM  sace  dp.  (3) 

o 

If  the  surface  of  the  earth  is  assumed  to  be  black,  (3)  can  be  written  in 
the  integral  form 

.A 

Kv>e)    =  B[v,T(po)]   exp[-  |^     kCv.p^qCp1)    sece  dp'] 

o 

Po  E 


+ 


I  C  B[v,T(p)]   exp[-  |^  kCv.p^qCp1)    sece  dp']  (A) 


x  k(v,p)q(p)    sece  dp   , 

where   the   subscript  0  refers   to  the    lower  boundary   (surface).      This   form  of 
the  equation  shows   the   two  components  of  the   radiances:      that  arising   from 
the   surface   and  attenuated   by  the   atmosphere,    and   that  arising  from  the 
atmosphere.      It  may  be  noted   that  a   "black"  cloud   in  the    field  affects   only 
the   pressure   of  the   boundary   level. 

At  this  point  one  must  consider  the   instrument's    field-of-view.      If 
R(cd, cp1)    in   (1)    is  zero  outside   some   small  angle    from  the   optical  axis,    then 
sec0  varies   only  slightly  over  the    field  and   its  value   at   the   center  of   the 
field  may  be  used  to  represent  the  entire    field.      That   is,    integration  over 
6  is  not,    in  practice,    necessary. 

It   is  convenient  to  simplify  the   notation  by  introducing  the    fractional 
transmittance   of  the  beam  between   level  p  and   the  effective   top  of   the 
atmosphere 

P 
t(v>P>  ft)   =  exp[-  I     V  k(v,p')q(p')   sece  dp'].  (5) 

o 

Then    (4)   becomes 

Kv,e)    =   B[v,T(pJ]    T(v,Pn,  9) 


s 


t(v,p0,  e) 

B[v,T(p)]   dT(v,P,e).  (6) 

1 


Equation   (6)   expresses  the   radiance   for  a  single   spectral   frequency.      The 
finite  width  of  the   spectral  interval  being  observed  requires  an  appropriate 
integration  over  the   response   of  the   instrument.      Considering  the   case  of  a 
filter-type   radiometer,    the  measurable   radiance    for  the   jtn  filter  is 


r«- 


Kvj,e)   =  J  Kv,  9)   tjCv)  dv 


B[v,T(po)]    t(v,P0,  G)    t    (v)   dv 


r°^  T(v,Prt,e) 


JJ  B[v,T(p)]   t    (v)   dT(v,p,  0)   dv  1  (7) 

°  1  J  ) 

-f-J  t.(v)   dv 
o 

In  eq  (7),  t(v)  is  the  filter  transmittance . 

A  similar  integral  can  be  written  for  a  spectrometer. 

Equation  (7)  is  hopelessly  complicated  if  one  wants  to  retrieve  B(v,T). 
However,  as  pointed  out  by  Elsasser  (1938)  if  the  spectral  interval  over 
which  tj(v)  is  nonzero  is  small,  then  B(V,T)  varies  little  and  is  nearly 
linear  in  the  interval.   It  can,  therefore,  be  replaced  by  its  value 
B(v4,T)  at  a  properly  defined  mean  frequency  v<»  and  factored  out  of  the 
integral  with  respect  to  v  after  the  order  of  integration  is  interchanged. 
It  becomes  necessary  to  define  the  mean  transmittance  and  its  derivative  by 

\  T(v,p,  0)  t,(v)  dv 
t(vj,P,  6)  = 


^  tjCv)  dv 


(8a) 


and 


o 

oO 


5  (dT/dp)(»»p,  e)  tj(v)  dv 
tr  (vi'p'e)  =  c^ 

dp  \     t  (v)  dv 


(8b) 


In  the  instruments  discussed  in  this  technical  note,  the  spectral  inter- 
vals are  narrow  (about  10  cm"*)  and  the  filter  functions  ti(v)  are  nearly 
symmetrical  about  a  vertical  line  through  the  centroid.   Therefore,  one  may 
set  vi  equal  to  the  mean  frequency  at  the  centroid  of  ti(v)  and  use  (8a)  and 
(8b)  to  rewrite  (6)  in  the  form 


KvjS)  =  B[Vj,T(po)]  T(Vj,Po,e) 


-S 


>t(vj,p0>  e) 


B[Vj,T(p)]  dT(Vj,P,0)  (6r) 


Equations  (6)  and  (6')  are  identical  in  form,  but  the  latter  contains  the 
approximation  B(v*,T),  introduces  a  mean  transmittance,  and  assumes  that 
v-  ■  v-.   The  step  in  going  from  one  equation  to  the  other  is  not  trivial, 
but  yields  physically  realizable  functions. 

B.   Spectral  Regions 

In  the  design  of  an  indirect  sounding  system  it  is  necessary  to  review 
the  absorption  bands  which  occur  in  that  part  of  the  spectrum  where  measur- 
able self -emission  occurs  in  the  earth's  atmosphere.   Measurements  in 
selected  "windows"  of  low  absorption  are  also  required  to  examine  surface 
temperatures  and  the  influence  of  clouds  upon  the  outgoing  radiation.   Prac- 
tical considerations  limit  measurements  of  the  spectrum  to  a  range  of  about 
3.5pm  to  about  20  mm. 

Absorptions  significant  for  indirect  temperature  soundings  in  the  atmos- 
phere result  from  oxygen,  carbon  dioxide,  water  vapor,  ozone,  nitrous  oxide, 
and  methane.   Table  1  lists  the  various  constituents,  their  mixtures  (by 
volume),  and  their  absorption  bands;  the  locations  of  several  important 
"windows"  also  are  indicated.   Radiances  in  bands  of  gases  with  fixed  mix- 
tures can  be  used  to  deduce  temperature  profiles;  for  gases  with  variable 
mixtures,  satellite  measurements  may  be  used  to  estimate  the  mixtures. 
Figure  1  (Valley  1965)  displays  the  approximate  strengths  and  shapes  of  the 
bands,  but  not  the  fine  structure  of  the  individual  lines. 

Table  1  shows  that  temperature  profiles  must  be  deduced  from  measurements 
in  one  or  more  of  three  spectral  regions: 

The  combined  CO2  and  N20  bands  in  the  2160-2393  cm"1  region; 
the  CO2  band  between  550  cm"l  and  760  cm"1;  and 
the  oxygen  complex  in  the  1.7-2.3  cm-1  region. 

Determination  of  the  mixtures  of  ozone  and  water  vapor  may  result  from 
measurements  (estimates  of  these  may  be  needed  for  temperature  determina- 
tions) in  the  following  ranges: 


For  water  vapor,  1200-2000  cm"1,  near  800  cm"1,  or  0.74-550  cm"1, 
for  ozone,  1000-1100  cm"1. 

Table  1. --Absorption  bands  and  windows  important  to  indirect  soundings 


v 


Mixture 
X  by  volume 


(cm"1)  (urn)      Species      (percent)     Application 

0.00005     fair 

033    |  temperature 


2,160-2,270 

4.4-4.6 

N20 

2,240-2,393 

4.18-4.5 

C02 

560-   770 

13  -  18 

co2 

1.7-2.3 

4300-5900 
(51-69  GHz) 

°2 

1,200-2,000 

5.0-8.3 

H20 

0.74-  800 

12.5-13,500 

(22.2  GHz) 

H20 

1,000-1,100 

9-10 

o3 

650-  750 

13.3-15.4 

°3 

2,393-2,870 

3.5-4.18 

window 

770-1,000 

10.0-13.0 

window 

0.9-1.4 

7,000-11,000 
(27-42  GHz) 

window 

<  0.65 

> 15, 000 
(19.5  GHz) 

window 

f0. 031-0 
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variable     j  H20  abundance 
(0-2)      J       and  cloud 


(variable 
0-0.0003 


emissivity 

O3  abundance 

or  transmittance 

Surface  and 
cloud  temperature 

j        Surface 
j  temperature 


There  are  two  practical  ranges  for  estimating  water  vapor  concentration. 
These  ranges  are  determined  partly  by  instrumental  capabilities;  they  are: 
near  the  0.74  cm"1  line  with  microwave  radiometers,  and  on  the  low  wavenum- 
ber  side  of  550  cm"1  using  infrared  techniques.   The  lower  limit  for  most 
purposes  is  about  300  cm"1.   Below  that,  absorption  by  water  vapor  is  so 
great  that  the  earth's  radiance  arises  from  near  the  tropopause.   This 
causes  insensitivity  to  the  radiance  of  the  water  vapor  mixture;  that  quan- 
tity, then,  becomes  largely  indeterminate. 

Two  minor  absorption  bands  deserve  mention  because  they  coincide  with 
more  important  bands.   In  the  region  650-750  cm"1,  the  701  cm~l  ozone  band 
imposes  a  variable  absorption  upon  the  nearly  constant  absorption  by  carbon 
dioxide.   Strictly  speaking  the  ozone  mixture  should  be  known  if  one  wishes 
to  infer  the  temperature  profile.   From  experience  with  the  SIRS  instrument 
it  has  been  found  that  a  rough  estimate  of  the  ozone  usually  suffices  for 
this  purpose.   In  a  similar  way,  the  methane  band  (not  shown  in  table  1) 
between  1200  cm"1  and  1400  cm"1  has  been  found  by  Hanel,  et  al.  (1971b)  to 
render  that  part  of  the  v2  band  of  water  vapor  to  be  useless  in  retrieving 


the  mixture  of  that  constituent  (the  reverse  is  also  true  for  the  methane 
mixture) . 

C.   Gas  Transmittances 

Equation  (61)  shows  that  the  measured  radiance  depends  on  the  atmospheric 
transmittance  and  on  the  Planck  function.   To  calculate  the  Planck  function, 
it  is  necessary  to  know  the  transmittance  in  most  temperature  retrieval 
methods. 

The  spectral  transmittance  of  an  atmospheric  gas,  i,  between  pressure 
level  p  and  the  satellite  is  defined,  as  in  (5),  by 

1  rP 
TiUp.O)  -  exp[-  I  V  kiCv'.p')  q±    (p1)  sece  dp']  .  (5«) 

o 

It  follows  that  if  N  constituents  absorb  at  the  same  frequency,  the  total 
spectral  transmittance  is 

N 

t(v,p,6)  -  J]~    TiUp.e)  .  (9) 

i=l 

For  a  satellite  instrument,  with  its  finite  spectral  bandpass,  a  new  trans- 
mittance must  be  defined.  This  incorporates  both  the  optical  properties  of 
the  atmosphere  and  the  characteristic  of  the  instrument's  spectral  response, 
as  in  (7)  and  (8).  In  the  case  of  a  dispersed  spectrum  with  the  detector  at 
the  focal  plane,  the  transmittances  must  be  defined  by 

©o  o° 

Ti(Vj,p,e)   =  ^  T.(v,P,9)  f(vj-v)  dvA  fj(v)dV  (10) 

o  o 

where   v.   is  the  abscissa  of  the   centroid  of  the   instrument's   "slit   function" 
f,    for  the   jt"  spectral   interval.      Similarly,    the   transmittances   for  a 
"straight  through"   filter  type   of   instrument  are   given,    from  Eq   (8a),    by 

•^(vyP, e)  -  \  T^p, e)  tj  (v)  dv/  \  t.(v)  dv,  (li) 

o  o 

where  v.  is  the  abscissa  of  the  centroid  of  the  filter  transmittance, t. 
J 

In  applying  (10)  and  (11)  to  the  atmosphere,  it  is  usually  assumed  that  (9) 
may  be  applied  in  the  form 

N 

t(Vj,p,  s)  =  77  ^(v  ,P,e)  .  (12) 

i=l 
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This  may  be  justified  only  on  practical  grounds;  it  assumes  that  there  is 
no  systematic  relation  between  the  absorptions  by  any  two  gases.   It  is 
required  wherever  gases  with  variable  mixtures  are  involved.  Tests  of  the 
validity  of  this  assumption  must  be  made,  and,  if  necessary,  modifications 
made  to  suit  each  spectral  interval. 

Transmittances  for  the  fixed  constituents  can  best  be  made  by  point-by- 
point  calculations  through  spectral  lines  for  a  series  of  pressure  levels. 
Drayson  (1966)  has  performed  these  calculations  for  CO2  using  line  positions 
and  parameters  taken  from  laboratory  experiments  and  from  theory,  for  a 
spectral  mesh  which  is  much  smaller  than  the  line  widths.   When  the  results 
were  used  with  SIRS,  a  slight  modification  in  some  spectral  intervals  was 
required  to  bring  the  observed  temperature  profiles  from  radiosondes  into 
agreement  with  the  measured  radiances  through  eq  (61).   Similar  calculations 
have  been  made  in  the  2000-2500  cm-1  region  for  CO2  (McClatchey  1965)  and  in 
the  4-6  mm  region  for  O2  (Meeks  and  Lilley  1963). 

Gases  whose  properties  vary  require  a  different  approach.   Because  the 
point-by-point  calculations  consume  so  much  time,  simpler  approaches  are 
required.   One  method  is  to  compute  t  from  eq  (10)  or  (11)  using  point-by- 
point  calculations  for  a  series  of  homogeneous  path  lengths,  varying  pressure 
and  temperature.   The  transmittances  may  then  be  fit  to  a  suitable  expansion 
of  a  function,  W,  in  three  variables 

T(Vj,T£,  pm,  un)  =  expt-WjO^,  pm,  un)]  ,  (13) 

where  u  =  optical  path  length,  and 
^,m,  n  =  indices  for  temperature,  pressure,  and  optical  path  length. 

Then  any  atmosphere  may  be  fit  incrementally  for  quadrature  purposes  by 

6T(Vj,T,p,u)  «  -exp[-W(T,p,u)]  6W(T,p,u)  .  (14) 

For  water  vapor,  Smith  (1969a)  found  that  computed  transmittances  could  be 
fit  well  to  a  series  expansion  of  the  form 

xi  yi  zi 
W(T,p,u)  =  exp[rC  X  Y  Z  x)    , 
i 

where  X,  Y,  and  Z  involve  logarithms  of  u,  p,  and  T,  and  the  coefficients, 
C-,  were  determined  from  a  least-square  solution.  The  quantities,  x^,  y., 
z-,  take  on  values  of  0,  1,  and  2  in  several  combinations. 

Ozone  is  not  so  well  defined  in  terms  of  the  line  parameters,  so  no  point- 
by-point  calculations  have  been  performed  in  the  infrared.   However,  the 
701  cm"1-  band  is  rather  weak,  absorbing  only  about  5  percent  for  an  average 
ozone  amount  in  the  atmosphere.   From  the  laboratory  measurements  of  McCaa 
and  Shaw  (1968),  a  good  estimate  of  atmospheric  ozone  transmittances  can  be 
made.   No  temperature  adjustments  are  justified. 
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Transmittances  by  water  vapor   in  the    infrared   have   recently  been    found   to 
contain  an  element  which  is  not  explained   by  classical  approaches.      In  the 
•"windows",    for  example,    Varanasi  et  al.    (1968)    found   that   there    is  a   con- 
tinuum whose   absorption  coefficient   is   proportional   to  the   partial  pressure 
of  water  vapor,    and  which  has   a    large   negative   exponential  dependence   upon 
temperature.      Bignell    (1970)    found   the   same   phenomenon  in  the  windows   and    in 
the   rotation  band   above   200   cm"1.      Thus,    to  explain   the   observed   atmospheric 
transmittances   of  water  vapor   in  the   continuum  it   is   necessary  to  multiply 
the   right-hand   side   of  eq   (13)    by 


T(vj,P,T)    =  exp[-Kj   V  p«    q2   FoP   sece  dp'] 


(15) 


where  Ki  includes  the  constant  part  of  the  absorption  coefficient  and  other 
constants.   It  has  been  found  that  n=5  yields  reasonably  good  results. 

The  explanation  by  Varanasi,  et  al.  (1968)  is  that  this  absorption  is 
caused  by  the  water  vapor  dimer,  (H20)2.  Observations  fit  the  chemical 
equilibrium  theory. 

Eq  (13)  applies  for  the  lines  in  the  spectrum;eq  (15)  applies  to  the  con- 
tinuum.  Therefore,  atmospheric  transmittance  of  water  vapor  calculated  from 
considerations  in  this  section  now  becomes  the  product  of  eq  (13)  and  (15), 
in  order  to  include  the  effects  of  both  the  line  and  continuum  components. 

Eq  (6')  defines  the  relationship  between  the  satellite  measured  radiance, 
I,  and  the  atmospheric  temperature  profile,  expressed  through  B[T(p)].   How- 
ever, the  radiances,  I,  also  depend  on  the  atmospheric  transmittances,  t. 
It  is  therefore  evident,  from  the  above  discussion  of  transmittance,  that 
difficult  problems  can  arise  in  the  calculation  of  T  from  measured  I's  with 
methods  based  on  eq  (61). 

3.   INVERSION  PROCEDURES  -  CLOUDLESS  CONDITIONS 

The  integral  form  of  the  radiative  transfer  equation  for  a  cloudless  atmos- 
phere, eq  (6'),  is  given  in  Section  2.   That  equation  can  be  modified  to 
yield  the  relation 

x 
o 

I(v)  =  B[v,T(x0)]t(v,x0)  -  V  B[v,T(x)]  d^x>  dx  (16) 

"       u    «J  dx 

o 

An  "inversion  procedure"  or  more  precisely  the  "inverse  problem"  of  radia- 
tive transfer,  is  a  method  for  finding  the  temperature  T  in  eq  (16),  as  a 
function  of  the  independent  variable  x,  from  satellite  observations  of  the 
radiance,  I,  at  select  mean  (or  central)  frequencies  v.   The  variable  x,  can 
be  any  single  valued  function  of  pressure.   The  transmission  functions,  t, 
and  the  weighting  functions,  dr/dx,  are  assumed  known  in  some  inversion 
methods.   However,  even  if  I  is  completely  known  as  a  function  of  v.  ecl  (16) 
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Figure  1 . --Comparison  of  the  near-infrared  solar  spectrum  with  laboratory 
spectra  of  various  atmospheric  gases  (from  Valley  1965). 
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generally  does  not  have  a  unique  solution  (Courant  and  Hilbert  1953,  pp. 
159-160).   In  practice  I(v)  is  measured  at  only  a  finite  number  of  frequen- 
cies, v»  and  when  errors  in  the  observations  are  considered,  the  solution  to 
the  inverse  problem  may  be  unstable  in  the  sense  that  small  errors  in  I(v) 
yield  disproportionately  large  errors  in  T(x)  (see  Phillips  (1962)  for  a  dis- 
cussion of  this) . 

In  view  of  these  difficulties  it  is  apparent  that  temperature  accuracies 
will  vary  from  one  solution  algorithm  to  another.   For  this  reason  three 
classes  of  inversion  procedures  will  be  discussed;  they  include  regression, 
inverse  matrix,  and  direct  methods.   Of  these,  in  a  cloudless  atmosphere, 
the  regression  method  need  not  deal  with  the  physical  eq  (16);  therefore  the 
physical  parameters  are  not  required.   However,  in  the  cloudy  atmosphere, 
the  physical  parameters  still  may  be  required  even  with  some  versions  of  the 
regression  method. 

A.   Modifications  to  Equation  (16) 

Before  considering  the  various  solutions  of  the  inverse  problem,  one  can 
do  three  things  to  make  eq  (16)  more  tractable.   First,  a  climatological 
mean  T(x)  (or  a  forecast  profile,  when  available)  can  be  subtracted  from 
T(x)  to  reduce  the  problem  to  that  of  solving  for  the  deviation  "h"  of  T 
from  T.  With  a  forecast  profile,  "h"  will  usually  be  smaller  than  with  a 
climatological-mean  profile. 

Second,  since  generally  more  than  85  percent  of  the  energy  in  the  ll^m 
"window"  comes  from  the  boundary  term  in  (16),  this  measurement  and  the 
first  approximation,  T,  can  be  used  in  (16)  to  determine  a  good  first  approx- 
imation of  the  surface  temperature  Ts(x0),  for  cloudless  conditions.   Conse- 
quently, the  boundary  terms  for  all  the  spectral  frequencies  can  be  calcula- 
ted and  subtracted  from  eq  (16).   This  puts  (16)  into  the  more  classical 
form  of  an  integral  equation  of  the  first  kind.   Separation  of  the  surface 
boundary  from  the  atmosphere  also  takes  care  of  the  situation  in  which  a 
discontinuity  exists  between  the  ground  temperature  and  the  shelter-air- 
temperature. 

Third,  eq  (16)  is  nonlinear;  it  is  necessary  to  remove  the  simultaneous 
dependence  of  B  on  v  anc*  T  under  the  integral.   Also,  B  depends  nonlinear ly 
on  T.   It  is  possible  to  form  a  linear  equation  in  which  the  Planck  function 
has  been  approximated  by  the  Taylor  expansion, 

B(V,T)  -  BUT)  =  h  dB<v»f>  (17) 

dT 

where 

h(x)  =  T(x)  -  f(x)  .  (18) 
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The  quantity  dB(v,T)/dT  can  be   calculated  at  every  pressure    level  in  ad- 
vance for  each  frequency  once  T  has  been  chosen.      Other  approximations  are 
given  in  Wark  and  Fleming   (1966),    Fleming  and  Smith  (1972),    Smith,    Woolf,    and 
Jacob  (1970),    and   Smith,    Woolf,    and  Fleming   (1972). 

If  these   three   steps  are   incorporated   into   (16),    it  can  be  written  as  a 
linear  Fredholm  equation  of  the   first   kind: 


x 

o 


(v)   =  ^     K(v'x)  h<x>dx  ' 


(19) 


where 

x 


r(v)   =  Kv)    -  B[v,Ts(x  )]t(v,x0)   +  J     B[v,T(x)  ]d](v  x>   dx  (20) 

"  o  dx 

v(     v\     -  dT<v>x>    .    dB(v,T) 
MV,  X)    = -— y 

dx  dT 

and  h(x)    is  defined   by   (18). 

Since   I(v)>    and  hence   r(v),    is  measured  at  a   finite  number,    M,    of   frequen- 
cies,   and   the   kernel,    K(v,  x),    is  generally  known  only  in  tabular   form,    eq 
(19)    is  approximated   and   solved  numerically.      Partition  the   interval    [0, x0] 
at  N  +  1  points  and    let  w:,    j   =   1,2,  ...,N  be   the   quadrature  weights  associ- 
ated with  Xj,    then   (19)   can  be   approximated  by 

N 

ri   =  E     aii  hi    >      i   =   1>2,...,M  , 
j  =  l  J 

or  in  matrix  form  by 

r  =  A  h,  (21) 


where 


a.j  =  WjK(Vi,Xj)AXj 


The  quantities,  a..,  are  kernel  functions,  that  is,  they  are  values  by 
which  the  temperature  deviations  must  be  multiplied,  or  weighted,  to  get  the 
radiance  deviations  in  eq  (21).   The  shape  of  a^  is  mainly  determined  by 
dT/dx.   An  illustration  of  dT/dx  as  a  function  of  pressure  is  shown  in  fig- 
ure 2.   Generally  the  weighting  functions  arising  in  temperature  sounding 
problems  with  downward  pointing  instruments  are  very  smooth  and  overlapping 
(fig.  2),  with  the  result  that  the  matrix  A  in  (21)  is  usually  ill-conditioned 
with  respect  to  matrix  inversion.   Hence,  a  direct  solution  of  the  linear 
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equation  (21)  is  impractical  (see  Phillips  1962,  and  Wark  and  Fleming  1966). 
Other  methods  of  solution,  such  as  those  discussed  next,  must  be  used. 
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Figure  2. --Derivative  of  transmittance  with  respect  to  xCp^p2/7  (from 
Smith,  Woolf,  and  Fleming  1972). 


B.   Regression  Method 

If  the  transmittance  properties  of  the  atmosphere  are  not  known,  or  if 
there  is  a  large  uncertainty  in  the  accuracy  of  the  transmittance  functions, 
the  regression  method  should  be  used.  It  can  however,  also  be  used  even  when 
the  transmittance  functions  are  known  fairly  well.   This  method  also  has 
applicability  when  calibration  difficulties  or  measurement  biases  arise.   It 
was  used  operationally  with  SIRS  A  on  Nimbus  3. 

To  implement  the  method,  one  needs  radiosonde  measurements  that  are  coin- 
cident in  space  and  time  with  satellite  measurements.   Since  the  radiosonde 
measurements  do  not  extend  sufficiently  high  into  the  atmosphere  they  must 
be  extended  by  the  use  of  rocketsonde  or  grenadesonde  measurements;  in  some 
applications  the  use  of  climatological  data  is  adequate. 

Suppose  L  pairs  of  simultaneous  measurements  of  radiosonde  temperature  and 
satellite  radiances  are  available.   Let  the  matrix 


H  =  [hjk]  ,   j  =  l,2,  ...,N;   k=l,2,...,L, 


(22) 


represent  the  array  of  L  temperature  profiles  with  values  given  at  N  pre- 
scribed pressure  levels,  and  suppose  the  corresponding  L  sets  of  satellite 
measurements  are  given  by  the  matrix 
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R 


[rik]  ,   i=l,2,...,M;   k=l,2,...,L.  (23) 


The  quantity  "h"  is  defined  by  (18)  and  r(v)  is  the_deviation  of  I(v)  from 
the  mean  of  all  L  satellite  radiance  measurements,  I(v). 

The  principle  of  the  regression  method  now  depends  on  the  following  rea- 
soning.  Eq  (16)  shows  that  each  radiance,  I(v)>  depends  on  the  temperatures 
at  many  levels  in  the  atmosphere.   Therefore  we  can  assume  that  the  tempera- 
ture at  any  one  level  influences  a  few  of  the  radiances  significantly.  More- 
over the  temperature  at  one  level  in  the  atmosphere  is  often  correlated  with 
the  temperature  at  other,  distant  levels.   Therefore,  the  temperature  at  a 
particular  level  might  be  statistically  related  to  all  the  radiances.   With 
this  expectation,  it  is  then  required  to  find  a  set  of  coefficients  by  which 
the  radiances  can  be  multiplied  to  give  the  "best"  temperature  estimate  in 
a  least-squares  sense. 

The  regression  problem  is  that  of  finding  the  N  X  M  matrix  C,  which  best 
fits  the  relation 

H  =  CR 

in  the   sense   of   least   squares.      In  other  words,    find   the  matrix  C  that 
minimizes   the  Euclidean  distance 

N,L  M  2 

5(C)   =  E       (h.k  -  EC     r     ) 
j,k=l     Jk     i=l  ^   lk 

=   tr  {(H-CR) (H-CR)Tj,  (24) 

where  "tr"  and  the  superscript  "T"  represent  the  trace  and  the  transpose, 
respectively. 

The  solution  is  found  by  differentiating  (24)  with  respect  to  the  elements 
of  the  matrix  C,  which  results  in  the  normal  equation 

dF      T       T 

jg  =  2HR  -  2CRR  =  0  , 

and   implies   that  the  matrix  of  regression  coefficients   is 

C  =  HrVrV1   ,  (25) 

where  the  superscript  -1  indicates  the  matrix  inverse.   To  avoid  the  possi- 
bility of  the  matrix  RRT  being  singular,  one  should  choose  L  much  larger 
than  M.   For  any  particular  set  of  M  radiance  values  £,  the  N  temperature 
deviations  ^h,  comprising  the  solution,  now  can  be  obtained  from  the  equation 
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h  ■  Cr  ,  (26) 

which  when  added  to  T  yields  the  desired  temperature  profile. 

An  advantage  of  the  linear  regression  method,  especially  in  the  cloudless 
case,  is  that  it  does  not  require  explicit  knowledge  of  the  weighting  func- 
tions.  However  a  linear  relation  is  often  assumed  to  exist  between  h  and  r, 
when  in  fact  this  is  not  true. 

The  regression  method  is  useful,  especially  when  certain  physical  quanti- 
ties, such  as  atmosphere  transmittance  are  unknown.   However,  if  the  weight- 
ing functions  and  calibrations  are  known,  the  retrieval  methods  described  in 
the  next  two  sections  are  appropriate. 

C.   Inverse  Matrix  Methods 

Full  Statistics 

These  methods  have  the  advantage  that  simultaneous  measurements  of  r  and  h 
are  not  required.   Instead  climatological  measurements  of  h  and  the  errors 
e(v)  of  measurement  in  the  corresponding  I(v)  can  be  used  to  estimate  the 
required  statistical  quantities.   This  will  be  discussed  after  eq  (30). 
Write  eq  (21)  as 

r   =  Ah  +/e  (27) 

where  e  are  the  errors  of  measurement. 

Now  apply  eq  (27)  to  all  of  the  corresponding  pairs  of  the  L  columns  of 
data  in  eq  (22)  and  (23),  put  them  into  arrays  in  the  appropriate  manner 
(Smith,  et  al.  1970)  to  obtain  the  matrix  equation 

R  =  AH  +  E  ,  (28) 

where  E  is  a  known  matrix  of  errors  corresponding  to  R.   Eq  (28)  essentially 
defines  the  error  matrix,  since  R,  A  and  H  are  assumed  given. 

Substitute  (28)  into  (25)  and  assume  that: 

(a)  The  errors  £  are  uncorrelated  with  Jh; 

(b)  the  errors  e-^  are  random  variables  with  zero  mean;  and 

(c)  the  quadrature  errors  are  negligible  with  respect  to  je.  It 
follows  that  the  matrix  products  HET  and  EH*  are  small  relative  to  any  of 
the  other  products,  and  that  in  (26) 

C  =  (HHT)  AT  [A(HHT)AT  +  EET]    .  (29) 
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Compare   this  C  with  the   one   in   (25)   and  note   that   this  one   incorporates  the 
correlations   of  temperature  between  pressure    levels   through  the   quantity 
(HHT),    while   the   one   in    (25)   does  not. 

When  the   right   side   of   (29)    is  appropriately  divided  and  multiplied  by  the 
scalar    (L-l),    and   it   is  applied   to   (26),    then   (26)   has  the   form 

£  -   ShA*    [AShAT  +  Se]"\   ,  (30) 

where  S^  and  Sfi  are,  by  (18)  and  assumption  (b)  above,  the  covariance 
matrices  of  h   and  je,  respectively.   Other  derivations  of  equation  (30)  are 
given  in  Foster  (1961),  Strand  and  Westwater  (1968a,  1968b),  and  Rodgers 
(1970). 

It  was  noted  after  eq  (28),  in  connection  with  the  derivation  of  (30), 
that  knowledge  of  both  R  and  H  are  necessary  to  obtain  E.   However,  in  prac- 
tice S^  can  be  estimated  from  climatological  collections  of  radiosonde  meas- 
urements and  Secan  be  estimated  from  the  satellite  onboard  calibrations;  this 
eliminates  the  need  for  simultaneous  cloudless  measurements  of  radiances  and 
temperatures  which  are  difficult  to  obtain  in  large  numbers.   But,  of  course, 
eq  (30)  now  assumes  a  good  knowledge  of  the  transmittances  which  are  con- 
tained in  A. 

Minimum  Information 

By  making  various  simplifying  assumptions  about  the  covariance  matrices  in 
(30),  it  is  possible  to  consider  many  special  cases.   Perhaps  the  most  impor- 
tant of  these  is  the  "minimum  information"  case  of  Foster  (1961).   It  assumes 
that  the  elements  of  Jh  are  uncorrelated,  the  elements  of  e^  are  uncorrelated, 
and  their  corresponding  variances  o*fi  and  <7|,  respectively,  are  constants, 
that  is, 

Sh  =  ah  XN  '    Se  =  ae  XM  ' 

where   I,   is  the  J  X  J   identity  matrix.      Under  these  assumptions,    equation 
(30)    reduces   to 

h  =  AT[AAT  +  YI   ]"  r  ,  (31) 


where 


Y  =  o*/ou      . 


'e   '"h 
Foster    (1961)    calls   Y  the   ratio  of  noise   to  signal  power.      A  direct  deriva- 
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tion  of  (31),  given  by  Smith,  Woolf,  and  Fleming  (1972),  shows  that  it  can  be 
interpreted  as  a  perturbation  technique  which  minimizes  the  deviation  of  the 
retrieved  profile  from  a  profile  forecasted  from  dynamical  considerations. 
This  is  the  method  used  operationally  with  SIRS  B,  Nimbus  4  data. 

If  the  matrix  identity 

T   T       -IT        -It 
A1  [AAX  +  YEm]    =  [AXA  +  YIN]   A1 

is  applied  to  (31),  then  the  solution  is  that  of  Twomey  (1963,  1965).   How- 
ever, all  the  methods  following  from  eq  (31)  can  be  categorized  as  the 
Tikhonov  (1963)  method  of  regularization. 

The  minimum  information  solution  is  the  first  of  those  examined  thus  far 
that  is  essentially  independent  of  a  priori  statistical  data,  except  for 
some  estimate  of  cTh»   *c  no  longer  contains  the  covariance  matrices;  there- 
fore, T  in  equations  (18)  and  (20)  now  can  be  any  initial  profile,  not  nece- 
ssarily the  mean  profile.   This  may  be  advantageous  because  one  can  use  a 
forecast  profile  which  is  generally  more  accurate  than  a  statistical  mean  for 
the  initial  profile.   Experience  suggests  that  the  accuracy  of  solutions  from 
the  algorithms  discussed  in  this  report  are  dependent  upon  temperatures 
selected  for  the  initial  profile. 

Iteration 

On  the  other  hand,  the  failure  to  use  statistical  information  can  have  an 
adverse  affect  on  the  solution,  but  some  compensation  for  an  inappropriate 
Y  can  be  made  through  iteration. 

This  is  accomplished  by  making  the  following  replacements  in  equations 
(17),  (18),  (20),  and  (21): 

a.  T^n+1)  for  T, 

b.  T<n>  for  T, 

c.  h}0*1'  for  h, 

d.  r<n>  for  r, 

e.  A(n)  for  A. 

When  these  changes  are  applied  to  equation  (31)  it  takes  the  form 

^<n+l)  m   A(n)T[A(n)A(n)  +  ^j-l^n)  #  (32) 

The  iterative  process  is  continued  until  an  nQ  is  found  for  which 

RMS(r(no))  ^RMS(e)  .  (33) 
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Other  iterative  methods  are  discussed  next. 

D.   Direct  Retrieval  Methods 

These  methods  do  not  depend  upon  a  priorj.  statistics  or  parameters,  nor  do 
they  involve  matrix  inversions. 

In  the  "regression"  and  the  "inverse  matrix"  methods  described  above,  the 
solutions  depend  on  finding  a  suitable  set  of  coefficients  by  which  to  multi- 
ply the  observed  radiances.   This  is  accomplished  by  inverting  the  integral 
operator  of  eq  (16). 

In  the  "direct"  retrieval  methods,  to  be  discussed  below,  eq  (16)  or  eq 
(21)  is  retained  in  its  integral  form,  without  inversion.   The  method  depends 
on  iterative  calculations  of  the  radiances  from  eq  (21),  or  eq  (16),  until 
the  residual  from  the  observed  radiance  is  smaller  than  a  predetermined 
"noise"  value.   Sometimes,  in  contrast  with  the  regression  or  inverse  matrix 
methods,  no  set  of  coefficients  which  multiply  the  radiances  is  derived  in 
advance  (see,  for  example,  Chahine  1968). 

Typically,  a  "first  guess"  is  selected.   Then,  from  eq  (16)  or  eq  (21), the 
radiances  are  calculated.   The  differences  between  the  observed  and  calcula- 
ted radiances  are  noted.   Then  the  "first  guess"  temperature  profile  is 
adjusted,  according  to  a  pre -determined  algorithm  involving  the  differences 
or  ratios  of  the  radiances.   The  adjustment  is  made  in  an  attempt  to  bring 
the  calculated  radiances,  from  eq  (16)  into  closer  agreement  with  the  ob- 
served radiances.   This  procedure  is  then  repeated  as  often  as  necessary  un- 
til convergence  occurs  (if  it  does  converge)  to  a  pre-set  "noise"  level  for 
the  radiances. 

In  order  to  keep  the  solution  algorithms  dimensionally  correct,  it  is  help- 
ful to  consider  solutions  of  Planck  function  instead  of  temperature  directly 
(Fleming  and  Smith  1972).  Let  v0  be  a  fixed  reference  frequently  which  is 
essentially  at  the  middle  of  the  range  of  frequencies  being  used;  let 


b(x)  =  B[v  ,T(x)]  -  B[Vo,T(x)]  , 


o 


and 


dT(v,x)    dB(v,T)/dT   ,  (34) 

KU'X)  "   dx     '  dB(v0,^/dT 

Once  b(x)  is  determined,  the  solution  T(x)  is  calculated  from  the  inverse 
Planck  function  of  b  +  B(v  ,  T)  . 

Perhaps  the  earliest  efforts  at  solving  integral  equations  of  the  first 
kind  by  direct  methods  were  by  Landweber  (1951).   His  interation  formula, 
when  applied  to  the  discrete  form  (21)  with  "b"  of  (34)  used  in  place  of  "h" 
and  with  (34)  used  for  "K",  is 
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b(n+1>=a)ATr(n)      .  (35) 


Where  the  superscripted  quantities  are  those  described  in  items  (a)  through 
(d)  in  the  sub-section  on  Iteration,  above  and  the  scalar  to  is  a  relaxation 
factor  which  assures  convergence,    provided  that 

0<U)<2///ATA//      . 

In  most  instances  the  rate  of  convergence  of  (35)  can  be  improved  by  repla- 
cing the  scalar  u)  by  a  matrix  D.   For  example,  a  slight  modification  of  the 
solution  suggested  by  Smith  (1970)  yields  the  result 

(n+D     T  (n)  ,  v 

b      =  DATr  (36) 


where  D  is  an  N  X  N  diagonal  matrix  with  diagonal  elements 
M 
*i=l 


/ 
Ijj  =  l/(  E  atj)   ,   j=l,  ...,N  . 


On  the  other  hand,  Conrath  and  Revah  (1972)  have  suggested  that  the  diago- 
nal elements  of  D  be 

N    M 


dji  =  l/(l    T   aikaiJ}   * 

JJ    '  k=l  i=l     J 

These  forms  of  (35)  converge  provided  that 


||  in  -  data||<i    . 

A  paper  by  Petryshyn  (1963)  discusses  formulas  (35)  and  (36)  as  well  as 
many  other  iterative  schemes. 

T  T 

Although  in  equations  (35)  and  (36)  A  was  presented  as  being  linear,  A  is 

in  fact  nonlinear.   This  is  evident  from  definition  (34)  of  the  kernel 

K(v,  x)  and  the  corresponding  numerical  representation  by  the  matrix  A  in  (21) 

in  which  it  is  nonlinear  because  T  changes  with  each  iteration.  (It  is 

assumed  here  that  t  is  independent  of  T.)  Consequently,  the  matrix  A  in  (35) 

and  (36)  should  carry  the  superscript  n  to  indicate  that  it  is  the  nth  iterate 

of  the  matrix  A. 
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One    can  perform  the    iteration  on  A  as   a  single-stage    iteration  at  the   same 

time   as   the   iteration  on ^r,    or  one   can  do  a   two-stage   iteration   in  which    the 
r   is   iterated   until  b  converges    for  a   fixed  A.      Then  A  is   iterated  again,    and 

the    iteration  on   r   proceeds   as  before    for  the   new   fixed  A,    and   so   forth.      In 

practice   the   single-stage   process  appears  to  be   adequate  when  convergence   is 
determined   by    (33). 

A  somewhat   different   iterative   approach  (Fleming,    1972)    is  presented  by 
the   scheme 


M 

S  q-a2  r<n) 
fn+1)   i=l  aiaijri 
b 


(37) 


where 


and 


J       M      (n) 

=  aiJ6iJ 
i=l      J 


/  N  2 


=  dB(Vi>Tjn))/dT 
ij     dB(Vo,TJn',)/dT 


The  definition 

,   v  (n) 

differs  from  that  given  for  (21)  and  does  not  depend  upon  T   .   Convergence 

of  (37),  as  defined  by  nQ  in  (33),  generally  is  more  rapid  than  for  (35)  or 
(36). 

All  of  the  direct  methods  discussed  thus  far  solve  the  nonlinear  eq  (16)  as 
the  limit  of  a  sequence  of  linear  equations.   However,  there  are  methods  that 
are  explicitly  nonlinear.   Perhaps  the  most  familiar  of  these  is  that  by 
Chahine  (1968,  1970,  1972). 

Chahine  assumes  that  the  Planck  function  B.(p.)  at  some  pre-selected  pres- 
sure levels,  p.,  can  be  related  to  the  upward  radiances  by 

*n+1rn.^  m  =  »nrT/«  m  I(vi)  (38) 


B.   [T<Pj>]  =  Bi[T(Pj)] 


inU) 
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con- 


I(v)  is  the  measured  radiance  and  In(v>i)  is  the  computed  value  of  I(vi) 
after  the  nth  iteration  of  the  vertical  temperature  profile.  When  In(vi) 
verges  sufficiently  close  to  the  observed  values  of  I(v^)  the  value  of 
Bn[T(pO]  is  taken  as  the  solution;  from  that  T(p.)  is  calculated  as  the 
value  of  T  at  the  pressure  level  p..   In  the  process,  the  temperatures  at 
levels  between  the  selected  pressure  levels  are  obtained  by  interpolation. 
This  method  is  nonlinear  and  does  not  depend  on  a^  priori  statistics,  nor  on 
a  set  of  pre -determined  coefficients;  the  quantities  by  which  the  observed 
I's  are  multiplied  changes  in  every  iteration.   Some  comparative  studies 
with  this  method  (Conrath  and  Revah  1972)  have  shown  that  its  solution  tends 
to  become  somewhat  unstable  when  the  number  of  measured  radiances  increases 
from  7  to  16. 

Smith  (1970)  has  modified  Chahine's  method.   Smith  assumes  that 
.n+1  „    .    n 


B 


(v^p)  -  B  (V.,P)+  [I(Vi)  -  I  (v^l   ,  (39) 


where  again  I(v)  is  the  measured  radiance.   But  now  Smith  allows  the  B's  to 
be  computed  at  all  pressure  levels.   A  separate  set  of  B's  will  be  generated 
for  each  v.   From  the  B's,  corresponding  temperatures  are  computed.   Smith 
then  proposes  that  the  temperature  at  any  pressure  level,  p,  be  given  by 


(n+1)       N  (n+1)  /  N 

T     (p)  =  E  T     (v.,P)AT(Vi,p)  /  S  AT(v,,P) 
i=l        x  I    i=l 


At  the  surface  AT(v.>p)  1S  replaced  by  t(v,P-)>  here,  p  is  the  surface 

X   s  s 

pressure. 

A  final  approach  that  should  be  mentioned  is  derived  by  Backus  and  Gilbert 
(1967,  1968,  1970).  Details  will  not  be  given  here  because  it  appears  that 
no  one  has  applied  the  method  to  the  temperature  sounding  problem.   This 
method  also  is  very  useful  from  a  diagnostic  point  of  view  in  that  it  yields 
explicit  procedures  for  determining  the  vertical  resolution  inherent  in  the 
weighting  functions  and  for  determining  "trade-off"  diagrams  between  resolv- 
ing power  and  measurement  errors  (Conrath  1971). 

E.   Conclusions 

In  comparing  temperature  retrieval  methods,  one  is  always  interested  in 
obtaining  the  most  accurate  solution;  but  this  often  depends  upon  the  kinds 
of  a  priori  data  available. 

For  example,  we  have  seen  that  without  explicit  knowledge  of  the  weighting 
functions  one  must  resort  to  regression  techniques.   However,  the  accuracy 
that  one  can  obtain  from  regression  methods  depends  strongly  on  the  nature  of 
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the  data  used  in  compiling  the  regression  matrices.   Consideration  must  be 
given  to  such  things  as  the  degree  of  correlation  that  exists  between  the 
quantities  measured  and  those  to  be  predicted,  and  the  degree  of  homogeneity 
of  the  data  versus  the  number  and  the  degree  of  independence  of  the  measure- 
ments. 

On  the  other  hand,  if  one  knows  the  atmospheric  transmittance  functions 
and  their  derivatives  well,  then  he  may  employ  one  of  the  inverse  matrix  or 
direct  retrieval  methods.   Of  these  methods,  the  full  statistical  method, 
eq  (30),  may  be  the  most  accurate  (Fleming  and  Smith  1972)  if  sufficient  a_ 
priori  information  is  supplied.   However,  experience  suggests  that  the 
accuracy  of  the  initial  profile  is  important  in  assessing  the  accuracy  of  the 
other  methods. 

Indeed,  there  is  a  point  at  which  the  accuracy  of  the  initial  solution  is 
sufficiently  good,  that  is,  r  is  sufficiently  small,  so  that  the  covariance 
matrices  in  eq  (30)  are  superfluous  provided  the  atmospheric  transmittances 
are  known.  At  this  point  the  minimum  information  case,  eq  (31)  or  (32),  is 
adequate.   But  in  this  case  the  direct  retrieval  methods  also  are  competitive, 

Another  consideration  is  the  question  of  instrumental  noise.  All  the 
methods  described  are  stable  in  the  sense  mentioned  on  the  first  page  of 
Section  3,  but  each  method  propagates  noise  into  the  solution  in  a  unique 
manner.   The  effect  of  noise  on  the  stability  and  "smoothness"  of  the  solu- 
tion may  be  different  for  different  methods.   The  full  statistical  method  is 
probably  the  one  most  capable  of  introducing  true  fine  structure,  but  noise 
can  introduce  fine  structure  that  is  not  real. 

In  an  operational  situation,  where  temperature  retrievals  must  be  made  in 
a  time  period  much  shorter  than  the  time  between  successive  fields-of-view, 
computer  running  time  is  of  paramount  importance.   Therefore,  because  of 
their  simplicity,  the  direct  retrieval  methods  that  have  fast  convergence 
rates  are  superior.  When  the  matrix  inversion  methods  are  used  with  itera- 
tion, the  inversions  must  be  made  repeatedly,  due  to  the  approximation,  eq 
(17),  which  is  incorporated  into  the  matrix  A. 

One  cannot  make  a  definitive  choice  of  method  without  first  considering 
the  many  trade-offs. 

Section  5  reviews  the  accuracy  obtained  operationally  by  the  regression 
method  with  SIRS  A  data  and  by  the  "minimum  information"  method  with  SIRS-B 
data.   However,  no  organized  effort  has  yet  been  made  with  real  data  to 
determine  which  methods  yield  the  most  accurate  results. 

4.   CLOUDS  AND  AEROSOLS 

We   have    so   far  considered   the   theory  and  retrieval  methods   for   the   cloud- 
less  atmosphere   only.      However,    the   atmosphere   is  rarely  completely  cloudless 
over  areas   in   the    fields   of  view  of  present-day  satellite   sounding  instru- 
ments. 
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Clouds  are   highly  variable    in  amount  and  characteristics.      They  drastic- 
ally modify  the   transmittance   and  emittance   of  the   atmosphere.      Therefore, 
when  the   satellite   observes   the   atmosphere    in  the   presence   of  clouds,    the 
observed  radiances  are   affected   not  only  by  the   atmospheric   temperatures, 
but  also  by  the   clouds.      It   is   necessary,    then,    to  eliminate   or   to  account 
for  the  effect  of  clouds  on  the   observed   radiances   to  retrieve   the   atmos- 
pheric  temperature    from  the   sounding  instrument  observations. 

Two  basic   approaches  have   been   suggested    for   dealing  with  clouds.      One 
utilizes  measurements   of  a   single    field-of-view;    the   other  employs   high  spa- 
tial resolution  measurements   of  multiple    fields-of-view.      In  one  version  of 
the   single    field-of-view  method,    an  attempt   is  made   to  estimate  what   the 
observed  radiances  would   have   been  had  there   been  no  clouds   in  the    field-of- 
view   (Smith,    Woolf,    and  Jacob   1970).      The   other  version   (Rodgers    1970)    treats 
the   satellite   observations  as   an  ensemble   of  cloud  and  clear  areas,    and  esti- 
mates  the    "most   probable"  temperature   profile  which  exists   in  the   observed, 
cloudy  area.      The  multiple    field-of-view  method    (Smith   1969c)    specifies   the 
clear  column  radiances   from  spatially  independent  radiances   observed  over  a 
non-overcast   regime. 

A.      The  Cloud  Problem 

Clouds  generally  exist  within  the    fields   of  view  of  a  sounding  radiometer. 
In  this  case 


Kv)    «  T)Icd(v)   +  (1   -   T))Iclr(v)  <40> 

where  I(v)  is  the  observation,  T\   is  the  fraction  of  the  field  covered  by 
clouds,  IC(j(v)  is  the  average  radiance  from  the  cloud  covered  portions  of 
the  field  of  view,  and  Icir(v)  is  the  average  radiance  from  the  cloud -free 
portions  of  the  field. 

Many  types  of  clouds,  such  as  fields  of  fair  weather  cumulus,  stratocumulus, 
and  cirrus,  tend  to  be  horizontally  stratified  and  irregularly  distributed  in 
the  horizontal  direction.   Cumulonimbus  clouds,  on  the  other  hand,  have  a 
wide  range  of  vertical  dimensions.   Depending  upon  the  field  of  view  of  the 
sensing  instrument  and  the  weather  situation,  several  different  types  of 
clouds  at  different  heights  may  contribute  to  a  given  radiance  observation. 
A  complete  parameterization  of  clouds  in  a  radiative  transfer  model  for 
deducing  temperature  profiles  would  introduce  more  variables  than  there  are 
independent  pieces  of  information  in  the  observations  from  present-day  instru- 
ments. 

For  the  simplest  case  of  a  single  layer  of  clouds 

Icd(v)  =  F(v,pc,0)  +  {rc(v)Id(v,pc,0)  +  ec(v)B[v,T(pc)] 

+  Tc(v)I(v,pc,ps)J  T(v,pc>0)  >  (4i) 
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where  p   is  the  pressure  of  the  cloud  layer,  and  rc,  g  ,  and  t  are  the  re- 
flectivity, emissivity,  and  transmissivity  of  the  clouds,  respectively.   The 
term  F(v,  Pc,0)  is  the  radiance  emitted  at  frequency  v  by  the  atmosphere  above 
the  cloud,  i.e.  between  the  cloud  pressure  pc  and  zero  pressure.   Id(v>Pc>0) 
is  the  downward  component  of  the  radiance  at  the  cloud  level,  B[v,  T(pc)  ]  is 
the  cloud  Planck  radiance,  and  I(v, PC,PS)  is  the  upward  component  of  the  radi- 
ance at  the  cloud  base,  i.e.  between  the  cloud  pressure  pc  and  the  surface 
pressure  ps.   In  writing  (41)  it  is  assumed  that  the  cloud  pressure  is  some 
effective  mean  value.   The  second  term  in  eq  (41)  is  the  radiance,  originat- 
ing at  the  cloud,  which  arrives  through  the  overlying  atmosphere  at  the 
satellite. 

At  intermediate  wavelengths  in  the  infrared  (e.g.  lS^m),  the  reflectivity 
of  most  clouds  is  near  zero.  Assume  rc  =  0,  so  that  Tc  =  l~e  ,  and  eq  (41) 
reduces  to  the  simple  form 

Icd(v)  =  ec(v)  Ib(v)  +  [l-ec(v)l  Icir(v)  ,  (42) 

where  Ih(v)  is  the  radiance  associated  with  the  opaque  'black'  cloud  condi- 
tion (i.e.,  ec  ■  1) .   Substituion  of  (42)  into  (40)  gives 

Kv)  =  <v(v)  Ib(v)  +  [l-(*(v>]  Iclr(v)  (43) 

where  ry(v)  is  the  effective  cloud  amount,  He  (v) .   After  making  the  appropri- 
ate substitutions  for  1^  and  Icir>  one  obtains 

Ps 

I(v)  =  B[v,T(ps)]  T(v,Ps,0)  -  ^  B[v,T(p)]  dT(v,p,0)  dp 

0  dp 


a(v)  {B[ 


v,T(ps)]  t(v,Ps,0) 


-  B[v,T(p  )]  t(v,P,,0) 

Ps 


c  -     -c 


_  C  B[v,T(p)]  dT(v,P>0)  dl  (44) 

Pc  dP 

The  sum  of  the  first  two  terms  on  the  right  is  equal  to  Iclr.   Even  for  this 
simple  case  of  a  non-reflecting  single  layer  of  clouds,  the  transfer  equation 
is  nonlinear  in  a/   a^d  p  and  cannot  be  solved  for  the  temperature  profile  by 
using  the  linear  inversion  methods  discussed  in  Section  3.   One  approach  to 
the  solution  of  this  nonlinear  cloud  problem  is  to  first  estimate  the  equiva- 
lent clear  column  radiances  from  the  measured  radiances,  and  then  to  solve 
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for  the  temperature  profile  from  the  equivalent  clear  column  radiances. 

Single  Field  of  View  approaches 

It  follows  from  (44)  that  the  equivalent  clear  column  radiance  for  any 
field  of  view  is 

Iclr(v)  =  Kv)  +  a(v)Q(v,Pc)  (45) 

where  Q(v>Pc)  is  the  term  enclosed  by  braces  in  eq  (44). 

It  is  tempting  to  determine  the  cloud  parameters  <y(v) >  and  p   from  auxili- 
ary data.   For  instance,  the  fraction  of  cloudiness,  7],  might  be  specified 
on  the  basis  of  high  resolution  visible  or  infrared  cloud  picture  data.   How- 
ever, this  fraction  will  not  be  equivalent  to  the  effective  cloud  amount  & 
at  the  frequency  v  unless  the  clouds  are  opaque  at  all  wavelengths,   e.g., 
thick  water  droplet  clouds.   Probably  a  more  severe  problem  is  the  difficulty 
in  registering  the  cloud  picture  data  with  respect  to  the  field  of  view  of 
the  sounding  radiometer.   A  small  degree  of  misregistration  could  lead  to 
significant  errors  of  71.   The  determination  of  the  cloud  pressure  altitude 
pc  from  auxiliary  data  is  even  more  difficult.   Even  if  the  cloud  top  pres- 
sure could  be  specified  (e.g.,  from  high  resolution  infrared  window  data) 
this  may  not  be  sufficient  for  assigning  an  effective  value  to  pc.   Most 
severe  is  the  fact  that  the  temperature  profile  below  the  cloud  must  be  known 
in  order  to  use  the  auxiliary  cloud  information  to  determine  the  correction 
<y(v)Q(v>  Pc) .   Nevertheless,  this  has  not  been  tried  yet,  so  it  is  not  known 
what  effect  such  a  procedure  would  have  on  the  accuracy  of  temperature  re- 
trievals. 

The  cloud  parameters  and  subsequent  cloud  corrections  to  the  radiances 
can  be  estimated  directly  from  the  set  of  radiance  observations  if  they  have 
a  common  field  of  view.   (Smith  1969b,  Smith,  Woolf,  and  Jacob  1970).   Since 
in  practice  the  "first  guess"  temperature  profile  is  now  estimated  before- 
hand, the  correction  is  evaluated  by  an  iterative  procedure  as  follows.   Us- 
ing measurements  in  two  spectral  intervals,  v,  and  vo>  whose  radiation  is 
sensitive  to  clouds,  one  can  eliminate  #(v)  from  (45)  after  assuming  the 
cloud  emissivity  is  constant  within  the  spectral  domain  of  observation.   The 
result  is 

[Iclr<V  -  Kv1)]Q(v2,Pc)  -  tlclr(v2)  -  Kv2)]Q(v1,Pc)  =  0.    (46) 

The  quantities  Icir(v)  and  Q(v,  pc)  are  calculated  from  an  initial  estimate  of 
the  temperature  profile  and  as  a  function  of  cloud  pressure,  and  pc  is  the 
root  of  the  nonlinear  equation  (46).   Once  p  is  determined,  the  effective 
cloud  amount  <y  is  determined  by  substitution  of  Q(v,  Pc)  into  (45)  for  spec- 
tral region  v     or  ^         These  cloud  parameters  then  are  used  in  equation  (45) 
to  obtain  Ic^r(v)  for  the  remaining  spectral  intervals.   These  equivalent 
clear  column  radiances  are  used  to  make  the  next  estimate  of  the  temperature 
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profile,  and  the  process  is  repeated  until  convergence  is  obtained.  Multi- 
level versions  of  this  cloud  correction  algorithm  have  been  derived  (Smith, 
Woolf  and  Fleming  1972). 

The  above  algorithm  has  been  applied  to  SIRS  and  IRIS  observations  from 
the  Nimbus  3  and  4  satellites.   (Smith,  Woolf,  and  Jacob  1970,  Smith,  Woolf 
and  Fleming  1972,  Conrath  and  Hanel  1970).   It  has  been  found  that  the  solu- 
tion for  the  cloud  parameters  and  temperature  profile  below  the  cloud  depends 
upon  the  initial  guess.   However,  this  method  can  be  used  to  reduce  the  in- 
fluence of  clouds  on  the  solution  profile  above  the  cloud  level. 

Another  approach  to  the  cloud  problem  is  presented  by  Rodgers  (1970).   The 
method  is  based  upon  a  maximum  probability  estimator  of  the  atmospheric 
state.   This  method  was  used  to  process  data  from  the  Selective  Chopper  Radi- 
ometer (SCR)  on  Nimbus  4  (Ellis,  et  al,  1970).   It  has  the  capability  of 
yielding  temperatures  below  cloud  levels  because  of  their  correlation  with 
the  temperatures  above.   Rodgers  (1970)  shows  that  the  accuracy  can  be 
improved  if  information  from  other  sources  (e.g.,  surface  temperature,  cloud 
picture  data,  and  forecast  profiles)  are  included  in  the  statistical  formu- 
lation. 

Multiple  Field  of  View  Approach 

There  are  two  major  difficulties  that  result  from  the  single  field  of  view 
approaches  outlined  above:   (1)  since  clouds  do  not  generally  transmit  much 
radiance  from  below,  the  radiance  observations  contain  relatively  little  in- 
formation about  the  temperature  of  the  air  below  the  cloud  level,  and  (2) 
different  temperature -cloud  structures  yield  the  same  radiances  so  that  the 
range  of  acceptable  solutions  for  the  temperature  profile  obtained  from  the 
measurements  becomes  broader.   Both  of  the  above  limitations  can  be  modified 
with  observations  obtained  from  a  spatially  scanning  radiometer  possessing 
sufficiently  high  spatial  resolution  (linear  ground  resolutions  better  than 
about  50  km).  With  such  measurements  it  is  possible  to  infer  clear-column 
radiances  from  the  fact  that  the  scale  of  horizontal  variability  of  clouds 
is  much  smaller  than  the  scale  of  variability  of  the  temperature  profile. 
A  numerical  algorithm  utilizing  this  property  for  determination  of  clear 
column  radiance  is  presented  by  Smith  (1968,  1969c). 

In  the  multiple  f ield-of-view  approach,  it  is  assumed  that,  in  two  geogra- 
phically independent  observations,  the  temperature  profile  and,  therefore, 
Iclr(v),  remains  the  same.   It  is  also  assumed  that  the  pressure  at  the  cloud 
top  (pc)  and  cloud  properties  (rc(v),  ec(v) »  anc*  T  (v))  remain  constant  for  a 
given  measurement  frequency,  and  only  the  fractional  cloud  cover,  T|,  varies. 
Then  from  (40) 

t    (  ,n  _  Ii<v)  -  N*I2(V) 

Iclr^v)  -  _± f ,  (47) 

(1  -  N*) 

where   Ii(v)    and  I2(v)   denote   the   observations  of  two  adjacent   geographically 
independent    fie lds-of -view  and  N*  is   the   ratio  of  the    fractional  cloud  covers 
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of  the  two  fields  of  view  (TU/TU).   The  quantity  N*  is  determined  from  ob- 
servations in  an  atmospheric  window  spectral  region.  Applying  (47)  to  a 
window  region  (v  »  w;  assume  Iclr  (w)  =  B[w,T(ps)]); 

Ii(w)  -  B  [w,T(p  )] 

N*  =  — ± — 

I2(w)  -  B  [w,T(ps)].  (48) 

The  surface  temperature,  T(ps),  is  known  from  other  satellite  observations 
(Smith,  Rao,  Koffler,  and  Curtis  1970).   If  measurements  from  the  sounding 
radiometer  are  available  in  two  independent  window  regions  (e.g.,  3.7p,m  and 
11.0y,m)  then  the  surface  temperature  can  be  obtained  directly  from  the  zero 
of  the  nonlinear  function 


[I1(w1)I2(w2)  -  I2(w1)I1(w2)]  +  B[w2,T(ps)][I2(w1)  -  I1(w1)] 

-  B[wlfT(p8)][I2(w2)  -  I1(w2)]  =  0  (49) 

which  results  from  the  elimination  of  N*  in  equation  (48).   In  practice  small 
corrections  for  molecular  absorption  must  be  made  to  the  window  observations 
before  using  eqs  (48)  and  (49). 

The  solution  described  above  for  the  clear  column  radiance  is  independent 
of  any  a  priori  knowledge  of  the  temperature  profile.   However,  the  solution 
is  valid  only  when  the  geographical  variation  of  the  observed  radiance  is 
due  solely  to  a  variation  of  cloud  amount.   A  variation  of  other  cloud  prop- 
erties or  ground  and  atmospheric  temperature  produces  erroneous  values  of  N* 
and  Ic^r(v).   Therefore,  the  fields  of  view  of  geographically  independent  ob- 
servations should  be  spatially  adjacent  so  that  the  observed  radiance  varia- 
tion will  tend  to  be  due  only  to  cloud  amount  variations.   In  addition,  this 
process  increases  experimental  error  by  comparison  with  cloud  free  areas; 
therefore  an  average,  determined  from  many  pairs  of  observations,  must  be 
taken  to  achieve  an  adequate  signal  to  noise  ratio  of  the  clear-column  radi- 
ances used  in  determining  the  temperature  profile.   High  spatial  resolution 
and  contiguous  scanning  ensure  geographical  closeness  of  the  observations, 
increase  the  probability  of  clear  fields  of  view,  and  produce  a  large  number 
of  independent  estimates  of  clear-air  radiance  for  a  given  geographical  area. 

The  above  method  for  calculating  Icir  was  initially  tested  using  aircraft 
radiometric  observations  (Smith  1971).   The  results  suggest  that  the  clear 
air  radiance  contribution  to  the  observations,  in  the  presence  of  a  partial 
cloud  cover,  can  be  determined  with  the  accuracy  needed  for  calculating 
temperature  profiles  down  to  the  earth's  surface.   It  is  expected  that  this 
method  will  be  used  to  process  observations  from  the  VTPR  (ITOS)  and  ITPR 
(Nimbus  5)  temperature  sounding  radiometers. 


30 


B.   The  Effects  of  Aerosols 

Another  source  of  error  in  temperature  profiles  derived  from  satellite 
measurements  is  the  effect  of  absorption  and  scattering  by  aerosols.   Ideally, 
the  influence  of  aerosols  can  be  accounted  for  in  the  temperature  inversion 
process  if  the  vertical  profile  of  aerosol  concentration  size  distribution, 
and  the  relevant  absorption  and  scattering  coefficients  are  known.   Until 
such  information  is  available  on  a  global  scale,  the  effects  of  aerosols 
must  be  taken  into  account  using  model  concentration  profiles,  or  else 
should  be  neglected  entirely. 

Recently,  Stowe  (1971)  investigated  absorption  and  scattering  by  aerosols 
of  radiation  in  various  spectral  regions  currently  used  for  atmospheric 
sounding.   He  found  that  the  scattering  by  the  aerosols  is  small,  but  that 
aerosols  may  be  treated  as  an  additional  absorbing  constituent  which  modi- 
fies the  atmospheric  transmission  function.   For  the  earth  and  atmosphere 
together,  the  aerosols  change  the  emerging  radiances  by  less  than  1  percent 
in  most  cases.   However,  aerosols  may  have  a  significant  effect  on  the  radi- 
ances emitted  by  them  from  the  lower  atmosphere  (e.g.,  at  750  cm"l).   Any 
aerosol-related  errors  which  might  occur  in  the  derived  temperature  profiles 
would  result  from  a  failure  to  take  account  of  the  modification  of  the 
"weighting  function"  in  the  lower  atmosphere  caused  by  the  aerosols. 

Marlatt,  et  al.  (1971)  recently  measured  atmospheric  aerosol  concentra- 
tions and  size  distributions  as  well  as  sea  surface  temperatures  from  the 
NASA  CV-990  aircraft.   Using  the  aerosol  concentrations  and  nearby  radio- 
sonde temperatures,  they  were  able  to  deduce  the  error  of  the  equivalent 
blackbody  temperature  in  the  window  region  due  to  neglect  of  the  observed 
aerosol  absorption.   For  very  hazy  conditions  with  cloudless  skies,  the 
error  was  only  of  the  order  of  IK,  in  agreement  with  Stowe' s  (1971)  calcula- 
tions. 

To  study  the  influence  of  aerosols  on  the  temperature  retrieval  problem, 
measurements  of  the  aerosol  concentrations  and  IR  radiances  for  a  variety  of 
spectral  regions  are  needed. 

5.   OPERATIONAL  EXPERIENCE  -  ERROR  ANALYSES 

A.   SIRS  A  Experience  -  Regression  Method 

Nimbus  3,  which  carried  SIRS  A,  was  launched  in  April  1969;  about  one 
month  later  the  observations  from  SIRS  A  were  put  into  operational  practice. 
The  regression  method  was  used.   The  method  of  deriving  the  regression  co- 
efficients was  summarized  by  Smith,  Woolf  and  Jacob  (1970)  as  follows: 

"1)   Samples  of  upwelling  radiance  observations  and  coincident  conven- 
tional temperature  and  geopotential  height  profile  observations 
are  obtained  .... 

"2)   Corrections  for  clouds,  high  terrain,  or  hot  terrain  are  calcula- 
ted and  applied  to  the  observed  radiances.   The  corrected  radi- 
ances are  then  converted  to  brightness  temperatures  using  the 
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Planck  equation. 

"3   The  brightness  temperatures  and  conventional  data  are  grouped  into 
the  various  latitude  zones.   For  each  zone,  two  sets  of  regression 
coefficients  are  computed  for  both  temperature  and  geopotential 
height.   One  set  is  based  on  six  channels  of  brightness  tempera- 
ture (that  is,  the  669.3,  677.8,  692.3,  699.3,  706.3  and  899.3  cm"1 
channels)  while  the  other  set  is  based  on  all  eight  channels. 

"4   Finally,  the  six-  and  eight-channel  regression  coefficients  and 
sample  mean  values  for  each  zone  are  stored  for  subsequent  inde- 
pendent calculations  of  temperature  and  geopotential  height  pro- 
files from  SIRS-observed  radiances." 

To  compute  the  regression  coefficients,  the  Northern  Hemisphere  was  divi- 
ded into  3  latitudinal  zones;  the  Southern  Hemisphere  was  divided  into  2 
zones.   Each  sample  period  in  the  Northern  Hemisphere  was  two  weeks;  in  the 
Southern  Hemisphere  it  was  six  weeks.   The  influence  of  clouds  was  removed 
in  various  ways  to  obtain  the  clear  column  radiances.  When  all  the  SIRS 
channels  were  operating  properly,  a  two-cloud  model  involving  two  of  the 
radiances  was  often  used.   When  some  of  the  channels  began  to  deteriorate, 
a  one-cloud  model  was  used  to  preserve  the  remaining  radiances  for  tempera- 
ture retrievals  after  the  clear  column  radiances  had  been  computed. 

The  regression  method  has  the  advantage  that  no  explicit  information  about 
the  absolute  value  of  the  SIRS  calibration  was  required.   It  was  only  neces- 
sary that  the  calibration  should  not  change  over  a  period  of  a  few  weeks;  the 
on-board  calibration  indicated  that  this  was  so.   No  explicit  information 
was  needed  about  the  C0£  or  other  transmittances.   With  other  methods  the  re- 
quirement for  transmittance  information  causes  difficulties.   However,  the 
sample  mean  was  used  as  "first  guess".   It  might  have  been  preferable  to 
have  used  a  better  first  guess;  the  12-hour  temperature  forecast  used  with 
SIRS  B  data  might  have  been  a  better  first  guess  than  the  sample  mean,  and 
might  have  yielded  better  temperature  retrievals. 

B.   SIRS  B  Experience  -  "Minimum  Information"  Method 

SIRS  B,  on  Nimbus  4,  was  launched  in  April  1970  and  put  into  operational 
use  in  June  1970.  With  SIRS  B  the  "minimum  information"  retrieval  method 
was  employed  operationally  (Smith,  Woolf,  and  Fleming  1972)  with  an  inverse 
matrix  method.   A  12-hour  temperature  forecast  produced  routinely  by  the 
National  Meteorological  Center  (NMC)  was  used  as  a  "first  guess".   This 
"minimum  information"  method  was  used  because,  after  experimentation,  it  was 
found  that  the  full  statistical  inverse  matrix  method  did  not  offer  suffic- 
ient improvement  and  consumed  considerably  more  computer  machine  time.   As 
applied  to  SIRS  B  data,  the  eight  measured  radiances  were  "corrected"  for 
cloudiness,  high  terrain,  and  hot  terrain,  as  they  had  been  in  SIRS  A  (see 
Smith,  Woolf,  and  Jacob  1970).   The  "corrected"  radiances  were  converted  to 
brightness  temperatures,  which  in  turn  were  converted  to  Planck  radiances  at 
700  cm"1  (the  approximate  central  frequency  of  the  15^m  observational  inter- 
val). Then,  by  assuming  that  the  climatological  temperature  variance  does 
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not  change  with  height,  an  expression  was  derived  for  the  Planck  function  at 
100  pressure  levels  in  the  vertical.   Upon  inverting  the  Planck  functions  the 
temperatures  were  obtained  at  the  corresponding  pressure  levels. 

C.   Error  Analyses 

There  has  been  essentially  a  full  year  of  operational  experience  with 
SIRS  A,  and  more  than  one  year  with  SIRS  B.   The  results  can  be  compared  with 
radiosonde  measurements  to  judge  the  accuracy  of  temperature  retrievals  under 
operational  conditions.   It  should  however,  be  noted  that  radiosonde  measure- 
ments, as  well  as  temperature  analyses  based  on  radiosondes,  have  errors  also. 

SIRS  A 

Smith,  Woolf  and  Jacob  (1970)  compared  the  SIRS  A  derived  temperatures 
with  objective  analyses  of  radiosonde  data.   The  verification  statistics 
were  compiled  from  data  over  Western  Europe  where  SIRS  radiances  were  ob- 
served within  3  hours  of  the  standard  0000GMT  and  1200GMT  radiosonde  release 
times. 

The  RMS  differences  between  SIRS-derived  and  radiosonde  temperatures  varied 
with  sky  condition  and  with  season.   When  the  skies  were  judged  to  be  clear, 
the  RMS  differences  were  generally  between  1.5  and  2K  for  June  and  September. 
In  January  they  were  mainly  between  2  and  3K.  When  high  clouds  were  judged 
to  be  present,  the  errors  increased.   The  RMS  errors  with  high  clouds,  were 
about  2K  in  June;  in  January  the  error  was  3.5K  near  the  500  mb  level,  but 
somewhat  less  at  other  levels.   To  summarize,  the  RMS  temperature  "errors" 
for  SIRS  A,  in  degrees  K  were: 

June -Sept.  1969  January  1970 

Clear  High  Clouds  Clear    High  Clouds 

1.5-2      2  2-3     3.5  (near  500  mb), 

somewhat  less  at 
other  levels 

In  part,  the  increased  error  in  January  1970  was  attributed  to  large  instru- 
mental noise  which  appeared  after  November  1969  in  the  714.3  cm"1  channel. 
The  temperature  "errors"  discussed  here  are  for  temperatures  at  specific 
pressure  levels  in  the  atmosphere.   If  temperatures  averaged  over  layers  in 
the  atmosphere  (or  thickness  values)  were  used,  the  "errors"  would  be  smaller. 

Another  type  of  comparison  was  made  by  Hayden  (1971).   Hayden  compared 
SIRS  A  thickness  values  with  corresponding  values  from  radiosondes  in  the 
vicinity.   Thickness  values  are  closely  related  to  the  mean  temperature  in  the 
layer.   From  August  to  October  1969,  the  errors  in  thickness  derived  from  SIRS 
was  comparable  to  errors  in  radiosonde  values,  except  that  SIRS  values  were 
somewhat  poorer  below  300  mbs  when  high  clouds  were  present.   However,  after 
November  1969,  tTie  SIRS  A  thicknesses  differed  markedly  from  the  radiosonde 
values  below  300  mb,  especially  when  high  clouds  were  present.  Tresumably 
this  was  also  related  to  the  loss  of  radiances  from  channel  714.3  cm"1  in 


January  1971 

All  Cases 

"Clear" 

High  Cloud 

2.4 

1.8 

2.5 

1.9 

1.2 

1.9 

1.6 

1.1 

1.6 
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SIRS  A  retrievals. 

SIRS  B 

Error  analyses  have   also  been  made   using  SIRS  B  operational  experience 
also.      Such  an  analysis,    for  January   1971,    appears   in  Smith,    Woolf,    and 
Fleming   (1972). 

The   SIRS  B  operational  retrievals  used   the    12-hour  NMC   forecast  as  a   first 
guess.      These   SIRS   retrievals  were    later  compared  with  the   NMC  analyses  which 
are   based  mainly  on  radiosonde   data.      The   resulting  RMS  errors,    in  degrees  K, 
were : 

July-Sept.  1970 

Pressure    (mb)  All  Cases 

200  2.5 

300  1.3 

500  1.1 

These  results  suggest  that  the  SIRS  B  temperature  profiles  may  have  been 
slightly  better  than  the  SIRS  A  results. 

The  analysis  for  July- September  1970  was  supplied  by  Brodrick  (manuscript 
in  preparation).   Brodrick  also  showed  that  if  the  NMC  analysis  itself  was 
used  as  a  first  fuess  the  "errors"  were  reduced  to  about  1.0K.   The  reason 
why  the  error  does  not  go  to  zero  in  that  case  is  presumably  that  the  NMC 
analysis  itself  does  not  give  the  vertical  temperatures  exactly  and  that  the 
SIRS  radiances  have  errors  in  them.   These  SIRS  radiances  contain  not  only 
random  instrumental  errors,  but  also  errors  introduced  by  clouds  and  other 
atmospheric  contaminants;  systematic  instrumental  errors  of  the  SIRS  may  also 
contribute  to  this  discrepancy. 

This  result  suggests  that  when  a  good  "first  guess"  is  obtained  the  results 
improve.   For  then  the  radiances  need  provide  only  small  adjustments  to  the 
"guessed"  temperature  profile  to  achieve  the  final  profile.   As  the  "first 
guess"  departs  further  from  the  true  profile,  larger  errors  may  be  found  in 
the  retrievals. 

6.   SOUNDING  INSTRUMENTS 

A.   Infrared 

Remote   sensing  of  atmospheric   temperature   profiles   from  spacecraft  was 
initiated  with  a   grating  spectrometer    (SIRS)   and  an   interferometer   (IRIS)   on 
the   Nimbus   3   spacecraft.      The  multi-detector,    fixed-grating,    Fastie-Ebert 
spectrometer,    an  early  version  of  the   SIRS,    was  designed   in   1960   (Dreyfus 
and  Hilleary   1962),    and  used  on  high  altitude   balloon   flights   in   1964 
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(Hilleary,  et  al.  1966)  and  in  1965  (Saiedy  and  Hilleary  1967)  to  test  the 
feasibility  of  obtaining  temperature  profiles  by  remote  sensing.   The  Fastie- 
Ebert  Spectrometer  optics  were  changed  slightly,  an  eighth  detector  was  added, 
and  its  electronics  were  redesigned  for  the  Nimbus  mission  (Wark,  et  al.  1969). 
A  grating  spectrometer  was  chosen  for  initial  sounding  attempts  because  of 
the  spectral  stability,  purity,  and  registration  which  it  affords.   The  use 
of  aluminum  for  both  the  Ebert  mirror  blank  and  the  supporting  framework,  and 
the  use  of  Pyrex  for  the  grating  substrate,  made  the  spectrometer's  wave- 
length calibration  and  spectral  bandpass  characteristics  essentially  indepen- 
dent of  instrument  temperature.   At  the  expense  of  a  large  field-of-view 
(12°),  no  fore-optics  were  used  between  the  radiometrically  balanced  beam 
modulator  and  the  earth  to  enhance  measurement  accuracy.   The  modulator 
referenced  the  measurements  against  the  radiance  of  intersteller  space  which 
is  essentially  zero  K. 

The  spectrometer  performed  well  in  space.   Its  calibrations  were  very 
stable  and  the  spectral  radiance  equivalents  of  the  RMS  signal  noises  ranged 
from  0.1  mw/(m2  sr  cm-1)  in  the  Hum  window  channel  to  0.3  mw/(m  sr  cm"1)  in 
the  15|jm  Q-branch  channel  in  spectral  intervals  5  cm"-*-  wide.   Limitation  of 
the  spectrometer  design  to  sensing  of  radiances  in  only  eight  discrete  spec- 
tral intervals  and  the  large  field  of  view  resulted  in  a  relatively  low 
volume  of  data  and  comparatively  simple  data  processing. 

After  about  two  months  in  orbit,  two  of  the  SIRS  thermistor  bolometer  detec- 
tors became  quite  noisy,  but  operational  use  of  the  data  continued  for  about 
one  year  until  similar  data  became  available  from  an  improved  SIRS  flown  on 
Nimbus  4. 

The  Nimbus  3  Michelson  interferometer  (Hanel,  et  al.  1970)  posed  severe 
engineering  problems.   Precise  control  of  its  temperature  was  necessary  to 
maintain  optimum  optical  alignment.   Prelaunch  testing  was  difficult  because 
of  its  sensitivity  to  vibration.  Many  significant  optical  engineering  prob- 
lems were  solved  during  an  intensive  and  sophisticated  development  program. 
Its  performance  in  orbit  was  compromised  somewhat  because  instrument  tempera- 
ture remained  about  5K  greater  than  expected,  but  a  wealth  of  earth  radiance 
spectral  measurements  (400  to  1600  cm"-*-)  were  obtained.   RMS  noise  equivalent 
spectral  radiances  between  0.5  and  1  mw/(m2  sr  cm""1)  were  achieved  with  5  cm"1 
resolution.   Its  field-of-view  was  8  degrees. 

The  Nimbus  4  spacecraft  carried  improved  versions  of  the  grating  spectrom- 
eter (Wark,  et  al.  1970),  and  the  Michelson  interferometer  (Hanel,  et  al. 
1971a),  and  also  a  Selective  Chopper  Radiometer  SCR,  (Abel,  et  al.  1970)  de- 
veloped in  the  United  Kingdom. 

The  Nimbus  4  spectrometer,  SIRS  B,  has  six  additional  thermistor  bolometers 
to  sense  spectral  radiances  in  the  18-  to  37|j,m  water  vapor  absorption  region, 
and  steps  its  sensing  beam  to  each  side  of  the  orbital  track  to  improve  the 
geographical  distribution  of  the  temperature  soundings.   Redesign  of  the 
optics  and  improved  thermistor  bolometer  detectors  permitted  the  SIRS  B  to  be 
smaller  and  lighter  than  the  Nimbus  3  spectrometer;  the  12  degrees  field-of- 
view  was  retained.   The  performance  of  the  CO2  band  and  ll-^m  window  sensing 
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channels  were  equivalent  to  those  of  the  Nimbus  3  instrument.   Noise  equiva- 
lent radiances  of  the  water  vapor  channels  ranged  from  0.25  to  0.45  mw/(n/ 
sr  cm"l).   Poor  electrical  connections  in  some  of  the  electronic  synchronous 
demodulator  circuits  produced  signal  range  saturations  and  spasmodic  gain 
changes,  but  operational  use  of  the  data  has  been  maintained,  except  for  the 
period  September  1  through  October  18,  1971,  and  continues  into  1972. 

The  Nimbus  4  Michelson  interferometer  achieved  a  spectral  resolution  of 
2.8  cm~l  over  the  spectral  range  400  to  1500  cm~l.   Its  spatial  resolution 
was  improved  to  94  km.   Its  noise  equivalent  spectral  radiances  are  about 
0.5  mw/(m^  sr  cm""l).   Problems  occurred  in  mounting  the  cesium  iodide  beam 
splitter  originally  proposed  and  a  potassium  bromide  component  similar  to 
that  used  in  the  Nimbus  3  instrument  was  substituted. 

The  selective  chopper  radiometer  uses  C0£  filled  cells  with  interference 
filters  for  spectral  selection.   It  is  able  to  provide  temperature  profiles 
extending  upward  to  50  kms,  although  calibration  problems  encountered  in 
orbit  may  affect  the  temperature  retrieval  accuracy. 

The  instruments  mentioned  above  can  be  considered  first  generation  atmos- 
pheric temperature  sounders.   They  employed  thermistor  bolometer  detectors 
and  achieved  effective  spectral  resolutions  of  about  5  cm"-'-  at  the  expense 
of  poor  spatial  resolution.   The  associated  satellite  experiments  proved  the 
feasibility  of  determining  atmospheric  temperature  profiles  from  satellite 
measurements.   However,  the  large  field  of  view  often  included  clouds,  so  the 
accuracy  of  temperature  retrievals  was  reduced.   Research  studies  indicated 
that  improved  soundings  might  be  obtained  in  partly  cloudy  regions  by  using 
statistical  techniques  if  the  instrument  field-of-view  can  be  reduced  and 
the  field  scanned  as  discussed  below.   It  was  also  realized  that  broader 
spectral  intervals  than  those  used  in  SIRS  could  be  used  in  the  CO2  15^™  band 
P-branch  and  R-branch.   A  multi-channel  temperature  profile  radiometer  (ITPR) 
(Smith,  Hilleary,  et  al.  1972),  and  a  sequential  filter  temperature  profile 
radiometer  (VTPR)  (Falbel  and  Zink,  1971)  have  been  specially  designed  for 
Nimbus  E  and  ITOS  spacecraft,  respectively,  to  exploit  those  ideas  and  to  pro- 
vide operational  weather  data.   Recent  developments  of  high  resolution  infra- 
red interference  filters  (Pidgeon  and  Smith  1964)  and  the  triglycine-sulfate 
pyroelectric  infrared  detector  (Putley  1970)  made  these  instruments  feasible. 
The  statistical  techniques  require  the  sensing  of  radiances  from  nearly 
adjacent  small  geographical  areas  (e.g.,  50  km  in  diameter)  preferably  arran- 
ged in  a  geographical  grid.   (see  Section  4,  Multiple  Field-of-View  Approach) 
The  measurements  of  radiance  from  adjacent  grid  elements  must  be  quite  inde- 
pendent. 

Both  the  ITPR  and  VTPR  scan  their  sensing  beams  perpendicular  to  their 
orbital  tracks.   The  scans  are  performed  as  a  series  of  steps  and  dwells, 
and  the  signals  are  sampled  and  integrated  during  the  scan  mirror  dwell  peri- 
ods.  On-board  calibration  tests  are  provided  for.   Certain  characteristics 
of  the  first  and  second  generation  instruments  are  listed  in  table  2. 
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The  Nimbus  E  instrument,  the  Infrared  Temperature  Profile  Radiometer  (ITPR) 
is  a  development  of  the  5-channel  medium  resolution  infrared  scanning  radi- 
ometer used  on  the  Nimbus  2  and  3  spacecraft.   A  preliminary  model  was  used 
on  an  aircraft  to  test  the  feasibility  of  deducing  "clear  air  radiances"  of 
a  partially  cloudy  atmosphere  (Smith,  Hilleary,  et  al.  1972).   The  ITPR  has 
been  completely  redesigned  for  use  on  the  Nimbus  E  spacecraft.   It  now  con- 
tains six  small  Cassegrainian  type  telescopes  which  share  a  common  45  degree 
scanning  mirror.   A  six  port  reflecting  rotary  chopper  near  the  six  field 
stops  permits  the  detectors  to  measure  the  differences  in  spectral  radiance 
between  the  earth  and  small  internal  blackbody  cavities.  A  dichroic  beam 
splitter  behind  one  telescope  reflects  the  near-IR  radiation  so  that  the 
3.8-um  and  lLl-ym  spectral  channel  detectors  can  share  one  telescope.   Inter- 
ference filters  are  used  to  select  the  spectral  intervals  shown  in  table  2, 
which  also  lists  the  noise  equivalent  radiances  achieved  during  tests.   Field 
lenses  are  used  to  image  the  Cassegrainian  secondary  mirrors  on  the  detector 
elements.   The  radiometric  characteristics  are  improved  (1)  by  controlling 
the  temperature  of  the  detector  housings,  filters,  and  chopping  reference 
blackbody  cavities,  (2)  by  telemetering  the  temperatures  of  the  scanning  and 
telescope  mirrors,  and  (3)  by  making  the  effective  detector  area  ring  shaped 
so  that  it  cannot  "see"  that  portion  of  the  telescope  aperture  obscured  by 
the  secondary  mirror. 

The  ITOS  Vertical  Temperature  Profile  Radiometer  (VTPR)  employs  a  single, 
larger  telescope  with  a  2.3  degree  field-of-view.   Its  eight  interference 
filters  are  mounted  in  a  wheel  behind  an  opaque  rotary  chopper.   The  single 
detector  views  first  the  earth  and  then  the  chopper  blade  through  each  filter 
in  turn  as  the  filter  wheel  rotates.   The  radiances  in  the  eight  spectral 
intervals  are  sensed  sequentially.   The  detector  housing  temperature  is  con- 
trolled and  the  temperature  of  a  shroud  surrounding  the  chopper  and  filter 
wheels  is  regulated  by  a  second  control  circuit  so  that  the  filter,  chopper, 
and  detector  temperatures  are  matched.   A  lens-and-cone  channel  condenser  are 
used  to  enhance  the  "light  grasp".   The  spectral  intervals  and  noise  equiva- 
lent radiances  of  the  VTPR  channels  are  listed  in  table  2. 

The  grating  spectrometer,  interferometer,  and  selective  chopper  design 
approaches  are  not  necessarily  superseded  by  these  new  designs.   Each  type 
of  instrument  offers  different  advantages  for  research  or  operational  pur- 
poses. More  advanced  spacecraft  grating  spectrometers  have  been  designed 
and  tested  on  balloon  flights  (Shaper  and  Shaw  1970).   Of  the  three  instru- 
ments described  the  selective  chopper  radiometer  (SCR)  approach  is  superior 
for  measurements  of  high  altitude  temperatures.   A  much  improved  16-channel 
SCR  is  being  developed  in  the  United  Kingdom  for  Nimbus  E.   An  improved 
interferometer  with  a  cesium  iodide  beam  splitter  functioning  satisfactorily 
is  aboard  the  Mariner  IX  spacecraft  in  orbit  around  Mars. 
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TABLE  2. --Infrared  Sounders 


Spacecraft  and 

Spectral 

interval: 

Spectral  interval 

Noise  equivalent* 

instrument 

center  frequency 
cm         <N/npj 

bandwidth 
cm"l 

spectral  radiance, 
mw/(m^  sr  cm"*-) 

Nimbus  3,  SIRS  A 

899.3 

11.2 

5.353 

0.11 

750.0 

13.3 

5.988 

.13 

714.3 

14.0 

5.35  7 

.20 

706.3 

14.2 

5.221 

.17 

699.3 

14.3 

5.102 

.17 

692.3 

14.4 

4.986 

.16 

677.8 

14.7 

4.748 

.21 

669.3 

15.0 

4.604 

.23 

Nimbus  3,  IRIS  B 

400  to 
1600 

25-6.2 

5 

0.5  to  1 

Nimbus  4,  SIRS  B 

899.0 

11.2 

7.36 

0.15 

750.0 

13.3 

6.83 

.19 

734.0 

13.6 

6.70 

.17 

709.0 

14.1 

6.13 

.21 

701.0 

14.3 

5.96 

.20 

692.0 

14.4 

5.81 

.16 

679.8 

14.7 

5.61 

.26 

668.7 

14.9 

5.42 

.20 

531.5 

18.8 

5.14 

.45 

436.5 

22.9 

5.00 

.41 

425.5 

23.5 

5.00 

.38 

302.0 

33.0 

5.00 

.29 

291.5 

34.3 

5.00 

.39 

280.0 

35.7 

5.00 

.42 

Nimbus  4,  IRIS  D 

400  to 
1500 

25  to  6 

.7    2.8 

.5 

ITOS,  VTPR 

835 

12.0 

8 

.25 

535 

18.7 

16 

.25 

747 

13.4 

10 

.25 

725 

13.8 

10 

.25 

708 

14.1 

10 

.25 

695 

14.4 

10 

.25 

677 

14.8 

10 

.25 

668.5 

15.2 

3.5 

.60 

Nimbus  E,  ITPR 

2675 
900 
745 
716 
691 
668.5 
505 

3.8 

11.1 
13.4 
13.9 
14.5 
15.0 
19.8 

450 
40 
19 
19 
19 
5 
80 

.003 

.19 

.19 

.19 

.20 

.60 

.20 

*Noise  Equivalent   Radiance    figures   calculated   from  pre-launch  test  data. 
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B.  Microwave  Radiometers 

A  microwave  temperature  profile  radiometer  is  planned  for  Nimbus  E  and 
another  for  Nimbus  F.   Both  radiometers  are  5-channel  instruments  with  one 
channel  in  the  22.235  GHz  water  vapor  band,  one  in  a  "window"  region  at  31.4 
GHz,  and  three  in  the  oxygen  complex.   The  performance  characteristics  are 
given  in  table  3. 

Table  3. --Nimbus  "E"  Radiometer  Performance  Characteristics 


H20 


Window 


0. 


0, 


Frequency  (GHz) 

Wavelength  (cm) 

RF  bandwidth  (MHz)  220 

Integration  time  (in  s-approx)    1.9 

ATrmS(°K) 

If   frequency  range    (MHz) 

Antenna  beamwidth  (deg)  10 


22.235        31.40 
1.349  0.955 


220 
1.9 
0.3  0.3 

10-110        10-110 
10 


53.65 
0.559 
220 
1.9 
0.7 
10-110 
10 


54.90 
0.546 
220 
1.9 
0.7 
10-110 
10 


58.80 
0.510 
220 
1.9 
0.7 
10-110 
10 


Antenna   types 


lens   loaded   horn 


Calibration  will  be   provided   by  two   loads,    one   at   the    instrument's   ambient 
temperature   and  one   that  will  be   radiatively  cooled   to  a   temperature   80C 
lower   than   instrument   ambient   temperature.      All  radiometers  will  be   of  the 
conventional  Dicke   type.      The   Nimbus   E   instrument  will  view  the   atmosphere 
directly  below  the    satellite   only.      The  Nimbus  F    instrument  will  be   a  very 
similar   instrument  except   that   is  will   scan.      The   scan  pattern  will  be  +  45 
degrees   from  nadir   in   13   steps   of  7.2  degrees  each.      The   frequencies   of  each 
of  the   02  channels  will  be  decreased    (i.e.,    shifted   to  a  slightly   less   absorb- 
ing  frequency)    to  compensate    for   the    longer   optical  path   lengths  encountered 
when  the   instrument   is   not  viewing  the   atmosphere   at  nadir. 

In  relation   to   infrared   instruments,    microwave   radiometers   offer  certain 
advantages   and   disadvantages.      One   very   important   advantage    is   their   greater 
cloud   penetration  capability    (Staelin   1969).      Another   is   the   ability  to  ob- 
tain weighting   functions   that   peak  in  the   stratosphere   and    lower  mesosphere 
(Lenoir,    1968).      This   can  be   accomplished  by  using  very  narrow  spectral   inter- 
vals  on   the   peaks   of   some   of  the  most   strongly  absorbing   lines. 

Some   disadvantages  can  be   attributed   to  the    fact   that   the   application  of 
microwave   radiometry  to   soundings    from   satellites   is   relatively  new.      In 
particular   calibration  techniques   are   not   as  extensively  developed   as   they 
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are  for  infrared  radiometers.   Thus  far,  microwave  radiometers  are  larger 
and  heavier  and  consume  more  power  than  infrared  radiometers,  particularly 
if  compared  on  a  per  channel  per  unit  field-of-view  basis.   However,  technol- 
ogy in  both  areas  is  developing  rapidly. 

Two  special  problem  areas  in  the  use  of  the  microwave  spectrum  are  variable 
emissivity  of  the  Earth's  surface  and  atmospheric  transmittance.   Surface 
emissivity  varies  over  a  considerable  range  of  values,  from  0.4  to  1.0 
(Staelin  1969),  and  can  be  affected  by  such  things  as  surface  roughness, 
surface  moisture,  etc.  Any  channel  whose  weighting  function  peaks  in  the 
lower  troposphere  will  be  affected  by  radiation  emitted  by  the  surface.  These 
factors  will  have  to  be  taken  into  account  in  experiment  design. 

Temperature  and  pressure  can  affect  the  atmospheric  transmittance  consid- 
erably. Most  inversion  schemes  require  an  accurate  knowledge  of  the  atmos- 
pheric transmittance.   To  supplement  the  earlier  work  of  Meeks  and  Lilly 
(1963),  Lenoir  (1968),  Lenoir,  et  al.  (1968),  and  Carter,  et  al  (1968),  a 
detailed  study  of  the  effect  of  temperature  and  pressure  broadening  on  O2 
absorption  is  being  carried  out  with  refraction  spectroscopy  (Thompson  and 
Vetter  1968,  Liebe  and  Dillon  1969).   Results  from  balloons  (Lenoir,  et  al. 
1968),  aircraft  (Carter,  et  al.  1968,  Rosenkranz,  et  al.  1971)  and  upward 
looking  radiometers  (Westwater  1972)  are  encouraging. 

C.   FURTHER  DISCUSSION 

The  optimum  instrumentation  system  for  atmospheric  temperature  soundings 
will  probably  include  infrared  and  microwave  radiometers.   The  spectral 
characteristics  of  optical  materials  and  detectors  will  not  permit  optimum 
design  of  a  single  radiometer  to  cover  the  entire  desired  infrared  spectrum. 
Exploitation  of  microwave  measurements  entails  a  different  engineering  disci- 
pline from  those  in  the  infrared  spectrum.   Several  radiometers  could  be 
designed  as  a  system  to  optimize  data  collection;  each  radiometer  in  the  sys- 
tem must  view  the  same  geographical  area  at  about  the  same  time.   The  radi- 
ometers should  be  calibrated  on  a  common  calibration  reference. 

The  absolute  accuracy  of  sounder  radiance  measurements  has  not  been  ade- 
quately emphasized.   All  of  the  sounders  discussed  here  have  included  means 
to  check  their  calibrations  in  orbit,  and  have  been  carefully  calibrated 
(Hilleary,  et  al.  1969)  in  vacuum-thermal  chambers.   But  comparisons  between 
SIRS  A  and  SIRS  B  grating  spectrometer  measurements,  and  uncertainties, 
(Hanel,  et  al.  1971b)  in  the  Nimbus  4  IRIS  calibration  indicate  significant 
problems.   Inter-comparisons  of  working  standards  and  mutual  consultations 
on  calibration  techniques  are  needed 

7.   SOME  REMAINING  PROBLEMS 

A.   Atmospheric  Transmittance 

Almost  all  retrieval  methods  depend  to  some  extent  on  a  knowledge  of  the 
atmospheric  transmittance  as  a  function  of  wavelength.   Most  operational 
experience  in  temperature  retrievals  involves  the  15-(jm  CO2  band.   As  was 
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pointed  out  in  Section  2,  there  was  a  systematic  difference  between  the  temp- 
eratures obtained  from  SIRS  A  radiances  and  the  radiosonde  temperatures.   It 
was  assumed  that  this  temperature  error  was  due  to  inaccuracies  in  the  values 
of  CO2  transmittance.   Therefore,  the  calculated  CO2  transmission  functions 
were  adjusted.   To  be  more  certain  about  the  CO2  transmittances,  more  measure- 
ments are  required  of  atmospheric  transmission  in  the  CO2  bands.   This  can  be 
done  by  observing  the  sun  in  the  15-^m  band;  the  atmospheric  emission  would 
then  be  a  minor  factor  by  comparison  with  the  hot  sun.   Suitable  platforms 
might  be  balloons  or  rockets;  some  progress  can  also  be  made  using  measure- 
ments from  mountain  tops. 

Such  a  mountain  top  experiment  has  been  performed  by  W.  L.  Smith  and  col- 
leagues.  They  mounted  an  ITPR  instrument  at  the  Mauna  Loa  Observatory  in 
Hawaii,  measured  the  energy  received  from  the  sun,  and,  separately,  the 
energy  received  from  the  adjacent  Earth's  atmosphere.   This  was  done  at  var- 
ious solar  elevation  angles,  for  wavelengths  in  the  CO2  band  at  15-ym  and  in 
the  rotational  water  vapor  band  near  20-|jm.   The  analysis  of  the  data  is  pro- 
ceeding. 

An  experiment  to  measure  the  atmospheric  transmittance  by  viewing  the  sun 
from  balloon  platforms  is  currently  being  planned  jointly  by  scientists  from 
the  United  States  and  from  the  United  Kingdom. 

Much  remains  to  be  done  on  the  problems  of  the  transmission  characteristics 
of  the  real  atmosphere  (mixture  of  gases).   For  the  temperature  retrieval 
problem,  it  is  necessary  to  estimate  the  concentration  of  the  variable  gases 
in  the  field  of  view.  Moreover,  since  the  transmittance  through  some  gases, 
such  as  water  vapor,  is  temperature  dependent,  a  non-linear  problem  exists; 
this  is  especially  important  for  the  lower  atmosphere.   To  aid  in  the  solu- 
tion of  this  problem,  the  temperature  dependence  of  the  transmittances  of  the 
atmospheric  gases  needs  to  be  investigated  further.   The  Mauna  Loa  experiment 
and  additional  laboratory  experiments  may  help  in  the  solution  of  the  problem. 

B.   Clouds  and  Particles 

Particles  may  influence  temperature  retrieval  for  the  lower  atmosphere 
(see  Section  4).   Some  numerical  experiments  with  simulated  "data"  are  being 
performed  (by  Wark,  Fleming  and  Stowe)  to  assess  the  magnitude  of  this  prob- 
lem. 

Although  transmittance  through  the  cloudless  atmosphere  still  poses  some 
problems,  clouds  remain  the  most  serious  source  of  errors  in  the  temperature 
retirevals. 

Until  now  actual  operational  experience  has  been  confined  to  instruments 
with  wide  fields  of  view  (SIRS  A  and  B) .   Section  4  discusses  methods  of 
handling  clouds  with  scanning  instruments  such  as  VTPR  and  ITPR,  which  will 
have  narrower  fields  of  view.   Measurements  from  these  instruments  should  be- 
come available  in  1972;  then  the  new  ideas  can  be  tried  out. 

However,  even  with  these  systems  further  problems  can  be  anticipated.   In 
some  areas  widespread  thin  cirrus  clouds  may  exist.   Since  microwave  measure- 
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merits  are  essentailly  unaffected  by  cirrus  clouds,  it  may  be  necessary  to 
use  microwave  measurements  to  cope  with  that  situation.   Microwave  measure- 
ments should  become  available  from  Nimbus  E  (1972)  and  Nimbus  F  (1974). 

The  "medium-resolution"  multiple-field-of-view  approach  of  VTPR  and  ITPR 
will,  it  is  expected,  yield  higher  overall  accuracy  than  the  SIRS,  wide- 
field-of-view  instruments.   However,  if  still  higher  accuracy  is  required,  it 
may  be  necessary  to  restrict  the  field  of  view  still  further,  and  to  attempt 
to  identify  those  areas  which  contain  no  clouds.   This  method  poses  instru- 
mental and  computational  problems.   Because  of  the  small  field  of  view,  very 
sensitive  detectors  would  be  required.   And  because  of  the  large  number  of 
observations,  expensive  data  handling  techniques  would  be  required.   It  is 
also  important  that  the  criteria  for  identifying  a  clear  column  are  suffici- 
ently definitive. 

For  these  methods  statistical  procedures  may  be  employed  to  select  the 
cloudless  areas  from  all  the  areas  viewed.   To  help  in  this,  it  may  be  use- 
ful to  observe  in  several  "windows"  such  as  the  infrared  windows  at  3.7ym 
and  ll|jm,  and  the  microwave  window  at  35  GHz.   However,  even  if  cloudless 
areas  cannot  be  identified  unambiguously,  the  technique  to  be  used  with  VTPR 
and  ITPR  can  be  employed;  many  more  sample  pairs  would  be  available,  and  the 
basic  assumptions  may  be  more  nearly  satisfied.   Therefore,  it  may  be  possible 
to  achieve  higher  accuracy. 

For  limited  areas,  the  geostationary  satellite  may  offer  some  advantages. 
Because  of  the  long  observational  dwell-times  attainable  from  a  geostationary 
satellite,  it  may  be  practical  to  seek  out  cloudless  areas,  and  observe  them 
for  long  time  intervals  to  reduce  instrumental  noise.  Moreover,  from  fre- 
quent observations,  it  may  be  possible  to  measure  the  time  variations  of  the 
vertical  temperature  profile;  this  may  be  important  in  certain  forecasting 
problems. 

C.   Retrieval  Methods 

In  Section  3,  various  temperature  retrieval  methods  were  discussed.   The 
regression  method  and  the  statistical  matrix-inversion  method  require  compu- 
tations of  "a  priori"  statistics.  A  fundamental  question  is  "how  should  the 
statistical  data  sample  be  selected?"  Specifically,  from  how  large  a  geo- 
graphical area  and  from  how  large  a  time  interval  should  the  samples  be 
drawn?  The  regression  method,  for  example,  is  sensitive  to  area  and  time 
changes. 

Some  of  the  methods  in  Section  3  derive  the  temperature  profiles  from  a 
form  such  as 


N 
T(p)  =  TQ(p)  +  £  c.(p)  Ali 

i 

in  which  the  c^'s  are  constants  derived  in  advance  from  statistical  data. 
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A  method  of  deriving  the  c. 's  might  be  the  following.   At  many  radiosonde 
stations  around  the  world  compute  the  c-^'s  for  each  station  from  past  statis- 
tical data.   Then,  when  satellite  observations  are  made  at  points  between  the 
radiosonde  stations,  interpolate  the  c^'s  to  the  satellite  position.   This 
might  work  better  than  large  area  statistics  applied  at  sub-satellite  points, 
especially  in  the  tropics  where  significant  temperature  changes  are  small. 

D.   Use  of  Other  Auxiliary  Information 

For  the  most  accurate  temperature  retrievals,  it  will  be  useful  to  employ 
all  available  information.   At  present,  for  example,  a  forecast  profile  is 
employed  to  give  a  "first  guess"  which  is  very  close  to  the  "true"  profile 
in  the  troposphere  and  lower  stratsophere.   In  addition,  satellite,  horizon- 
sensing,  radiance  observations  will  become  available  for  the  stratosphere 
from  Nimbus  F  in  1974.   This  may  yield  more  accurate  temperature  profiles  in 
the  stratosphere.   Preliminary  calculations  show  that  this  may  then  lead  to 
more  accurate  temperature  retrievals  in  the  troposphere  as  well. 

Cloud  information  can  be  obtained  from  auxiliary  satellite  observations. 
High  resolution  visible  and  infrared  "pictures"  may  help  in  specifying  cloud 
characteristics.   Such  cloud  information  may  be  useful  for  some  temperature 
retrieval  methods. 

Geostationary  satellite  cloud  motion  pictures  can  help  to  delineate  those 
areas  where  cirrus  clouds  exist.   Such  areas  might  be  handled  differently 
from  areas  which  do  not  contain  cirrus  clouds. 

In  some  areas,  clouds  and  cloud  types  are  related  to  the  temperature  struc- 
ture. If  correlations  between  such  clouds  and  temperatures  could  be  obtained, 
this  might  be  employed  in  Rodgers '  (1970)  method  of  temperature  retrievals. 

E.   Error  Analysis 

The  purpose  of  satellite  temperature  sounders  is  to  obtain  "accurate"  ver- 
tical temperature  profiles  over  the  whole  Earth.   How  should  "accuracy"  be 
determined?  The  verification  method  used  now  is  based  on  comparison  with  an 
independent  measurement  of  temperature.   The  only  "standard"  available  depends 
on  radiosonde  observations  for  the  troposphere,  and  radiosonde  and  rocket- 
sondes  for  the  stratosphere.   But  both  of  those  types  of  measurements  contain 
substantial  errors.  Are  analyzed  maps  better  "standards"  than  the  radiosondes 
themselves? 

In  utilizing  radiosonde  comparisons,  it  should  be  noted  that  the  precise 
location  of  the  satellite  field-of-view  is  probably  not  known  to  better  than 
30  km,  and  that  the  radiosonde  position  is  also  not  taken  into  account  pre- 
cisely.  In  winter,  when  large  horizontal  temperature  gradients  exist,  loca- 
tion uncertainty  is  a  source  of  error.   Moreover,  the  radiosonde  samples  a 
very  small  volume  which  contains  random  temperature  fluctuations,  while  the 
satellite  sounder  smooths  in  3  dimensions. 
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The  retrieval  temperature  accuracy  may  depend  on  several  factors.   For 
example,  what  is  the  influence  of  the  "first  guess"?  Experience  so  far  sug- 
gests that,  in  areas  with  good  radiosonde  coverage,  the  routine  temperature 
forecast  used  as  a  "first  guess"  produces  more  accurate  temperature  retrievals 
than  methods  which  use  climatology  as  a  first  guess.  More  work  needs  to  be 
done  to  verify  this. 

Error  estimates  have  appeared  in  many  theoretical  studies.  Among  the  fac- 
tors which  affect  "accuracy"  are  the  following: 

1.  Spectral  intervals  selected  in  the  4(jm,  15ypi,  and  microwave  bands, 

2.  instrumental  noise, 

3.  other  "noise"  (for  example  clouds),  and 

4.  uncertainties  in  estimates  of  atmospheric  transmittance. 

It  would  be  interesting  to  study  the  following  question:   "What  is  the 
highest  accuracy  attainable  from  the  use  of  all  known  techniques  -  instru- 
mental, observational,  and  mathematical?" 

Finally,  soundings  are  needed  primarily  to  improve  the  forecast  rather  than 
to  improve  analyses.   Experiments  are  needed  to  assess  the  impact  of  satel- 
lite soundings  in  forecasting  accuracy. 

F.   Sounding  Instruments 

Instrument  calibration  remains  a  problem.   The  absolute  calibration  in  the 
laboratory  involves  some  error  in  wavelength  determination  and  in  absolute 
radiance  at  the  estimated  wavelength.   These  two  factors  are  intertwined 
through  the  Planck  relation  for  blackbodies.   This  absolute  error  impacts 
temperature  retrieval  directly,  especially  in  non-regression  methods.   It 
will  impact  comparisons  especially  in  retrievals  between  two  different  instru- 
ments.  For  example,  SIRS  A  and  SIRS  B  observing  the  same  scene  do  not  neces- 
sarily yield  the  same  retrievals;  one  factor  is  the  calibration  uncertainty. 

Intercomparisons  of  calibration  would  be  even  more  difficult  between  infra- 
red and  microwave  instruments  since  these  involve  basically  different  types 
of  instruments  and  calibration  techniques. 

Instrumental  noise  inevitably  reduces  the  temperature  retrieval  accuracy. 
Work  to  reduce  the  instrumental  noise  by  various  methods  such  as  improved 
detector  sensitivity  will  doubtless  help  to  improve  retrievals. 
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